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Abstract. Resonance and decay phenomena are ubiquitous in the quantum 
world. To understand them in their complexity it is useful to study solvable 
models in a wide sense, that is, systems which can be treated by analytical 
means. The present review offers a survey of such models starting the classical 
Friedrichs result and carrying further to recent developments in the theory of 
quantum graphs. Our attention concentrates on dynamical mechanism under- 
lying resonance effects and at time evolution of the related unstable systems. 



Any general physical theory deals not only with objects as they are but also 
has to ask how they emerge and disappear in the time evolution and what one 
can learn from their temporary existence. In the quantum realm such processes 
are even more important than in classical physics. With few notable exceptions 
the elementary particles are unstable and also among nuclei, atoms and molecules 
unstable systems widely outnumber stable ones, even if the lack of permanence is a 
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relative notion — it is enough to recall that the observed lifetime scale of particles 
and nuclei ranges from femtoseconds to geological times. 

It is natural that the quantum theory had to deal with such temporarily existent 
objects already in its nascent period, and it did it using simple means suggested by 
the intuition of the founding fathers. As time went, of course, a need appeared for 
a better understanding of these phenomena even if there was no substantial doubt 
about their mechanism; one can cite, e.g., a critical discussion of the textbook 
derivation of the 'Fermi golden rule' |Fei lecture 23] in |RS[ notes to Sec. XII. 6], 
or the necessarily non-exponential character of decay laws |Ex[ notes to Sec. 1.3] 
which surprisingly keeps to puzzle some people to this day. 

The quest for mathematically consistent description of resonance effects brought 
many results. It is worth to mention that some of them were rather practical. 
Maybe the best example of the last claim is the method of determining resonance 
poles using the so-called complex scaling. It has distinctively mathematical roots, in 
particular, in the papers |AC71L IBC711 rSi79j . however, its applications in molec- 
ular physics were so successful that people in this area refer typically to secondary 
sources such as jMoQSj instead giving credit to the original authors. 

Description of resonances and unstable system dynamics is a rich subject with 
many aspects. To grasp them in their full complexity it is useful to develop a variety 
of tools among which an important place belongs to various solvable models of such 
systems. Those are the main topic of the present review paper which summarizes 
results obtained in this area over a long time period by various people including 
the author and his collaborators. As a caveat, however, one has say also that the 
subject has so many aspects that a review like this one cannot cover all of them; 
our ambition is to give just a reasonably complete picture. We also remain for 
the moment vague about what the adjective 'solvable' could mean in the present 
context; we will return to this question in Section [5] below. 

2. Preliminaries 

Before starting the review it is useful to recall some notions we will need fre- 
quently in the following. Let us start with resonances. While from the physics 
point of view we usually have in mind a single phenomenon when speaking of a 
resonance in a quantum system, mathematically it may refer to different concepts. 

We will describe two most important definitions starting from that of a resol- 
vent resonance. A conservative quantum system is characterized by a family of 
observables represented by self-adjoint operators on an appropriate state Hilbert 
space. A prominent role among them is played by its Hamiltonian H, or operator 
of total energy. As a self-adjoint operator it has the spectrum which is a subset of 
the real line while the rest of the complex plane belongs to its resolvent set g{H) 
and the resolvent z i— > {H — z)~^ is an analytic function on it having thus no sin- 
gularities. It may happen, however, that it has an analytic continuation, typically 
across the cut given by the continuous spectrum oi H — one usually speaks in this 
connection about another sheet of the 'energy surface' — and that this continuation 
is meromorhic having pole singularities which we identify with resonances. 

An alternative concept is to associate resonances with scattering. Given a pair 
{H,Ho) of self-adjoint operators regarded as the full and free Hamiltonian of the 
system we can construct scattering theory in the standard way ^AJSi, iRSj . in 
particular, we can check existence of the scattering operator and demonstrate that 
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it can be written in the form of a direct integral, the corresponding fiber operators 
being called on-shell scattering matrices. The latter can be extended to meromorhic 
function and resonances are identified in this case with their poles. 

While the resonances defined in the two above described ways often coincide, 
especially in the situations when H = —A + \^ is a Schrodinger operator and Hq — 
—A its free counterpart, there is no a priori reason why it should be always true; it 
is enough to realize that resolvent resonances characterize a single operator while 
the scattering ones are given by a pair of them. Establishing equivalence between 
the two notions is usually one of the first tasks when investigating resonances. 

In order to explain how resonances are related to temporarily existing objects 
we have to recall basic facts about unstable quantum systems. To describe such a 
system we must not regard it as isolated, rather as a part of a larger system including 
its decay products. We associated with the latter a state space H on which unitary 
evolution operator U : U{t) = e~*-^* related to a self-adjoint Hamiltonian H acts. 
The unstable system corresponds to a proper subspace T-L^ C % associated with 
a projection E-^^. To get a nontrivial model we assume that is not invariant 
w.r.t. U{t) for any t > 0; in that case we have ||£'u[/(t)^|| < for ifj £ Hu and 
the state which is at the initial instant t — represented by the vector ip evolves 
into a superposition containing a component in "H^ describing the decay products. 
Evolution of the unstable system alone is determined by the reduced propagator 



which is a contraction satisfying V{t)* — V{~t) for any t G M, strongly continuous 
with respect to the time variable. For a unit vector ^ e "Hu the decay law 



is a continuous function such that < Pjj){t) < Pip{0) = 1 meaning the probability 
that the system undisturbed by measurement with be found undecayed at time t. 

Under our assumptions the reduced evolution cannot be a group, however, it 
is not excluded that it has the semigroup property, V{s)V{t) = V{s + t) for all 
s,i e M. As a example consider the situation where "Hu is one-dimensional being 
spanned by a unit vector ip € H and the reduced propagator is a multiplication by 



where = Eh {—oo, A] is the spectral projection of H. If ip and H are such that 
the measure has Breit-Wigner shape, d{'ijj, E^ip) = ^ [(A — Ao)^ -I- |r^] ^dA for 
some Aq G M and F > 0, we get v{t) = e-»^ot-r|t|/2 giving exponential decay law. 
Note that the indicated choice of the measure requires a{H) — M; this conclusion 
is not restricted to the one-dimensional case but it holds generally. 

Theorem 2.1. |Si72| Under the stated assumptions the reduced propagator can 
have the semigroup property only if a(H) — M. 

At a glance, this seems to be a problem since the exponential character of the 
decay laws conforms with experimental evidence in most cases, and at the same 
time Hamiltonians are usually supposed to be below bounded. However, such a 
spectral restriction excludes only the exact validity of semigroup reduced evolution 
allowing it to be an approximation, possibly a rather good one. To understand 
better its nature, let us express the reduced evolution by means of the reduced 



V : Vit) = E^U{t) \ n 



■U 7 



(2.1) 



P^, : P^t) = \\V{t)^^ - \\EuU{t)^\ 
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resolvent, R^{z) :— E^Rh{z) \ Hu- Using the spectral decomposition of H we can 
write the reduced propagator as Fourier image, 

(2.2) V{t)^p= [ e-'^'dFx^p 

Jr 

for any ijj € 'Hu, of the operator- valued measure on M determined by the relation 
F{-oo, A] := E^E\"'\nn ■ By Stone formula, we can express the measure as 

the formula simplifies if the spectrum is purely absolutely continuous and the left- 
hand side can be simply written as -F'(Ai, A2). The support of F{-)ip is obviously con- 
tained in the spectrum of H and the same is true for suppi^ — Uv'e'H ^upp i<"(-)'0, 
if fact, the latter coincides with (j{H) 'Ex76' . 

In view of that the reduced resolvent makes no sense at the points ^ € supp F 
but the limits s-lim^_>.o+ R^i^ ± i'rj) may exist; if they are bounded on the interval 
(Ai,A2) we may interchange the limit with the integral. Furthermore, since the 
resolvent is analytic in p{H) = C \ (j{H) the same is true for i?^(-)- At the points 
of suppi^ it has a singularity but it may have an analytic continuation across it; 
the situation is particularly interesting when this continuation has a meromorphic 
structure, i.e. isolated poles in the lower halfplane. For the sake of simplicity 
consider again the situation with dimHu = 1 when the reduced resolvent acts as a 
multiplication by r'^iz) and suppose its continuation has a single pole, 

(2.3) rliz)^^ + f^,) 

Zp z 

for Imz > 0, where A ^ 0, f is holomorphic, and Zp Ap— iSp is a point in the 
lower halfplane. Since r^{X — iri) = r^{\+i'q), the measure in question is 

and evaluating the reduced propagator using the residue theorem we get 

(2.4) w(i) = Ae-'^p*-*pl*l - / e-*^*Im/(A)dA, 

^ Jr 

which is close to a semigroup, giving an approximately exponential decay law with 
F = 25p, if the second term is small and A does not differ much from one. At the 
same time, the presence of the pole in the analytic continuation provides a link to 
the concept of (resolvent) resonance quoted above. 

The main question in investigation of resonances and decays is to analyze how 
such singularities can arise from the dynamics of the systems involved. A discussion 
of this question in a variety of models will be our main topic in the following sections. 

3. A progenitor: Friedrichs model 

We start with the mother of all resonance models for which we are indebted 
to Kurt O. Friedrichs who formulated it in his seminal paper |Fr48] . This is not 
to say it was recognized as seminal immediately, quite the contrary. Only after 
T.D. Lee six years later came with a caricature model of decay in quantum field 
theory, it was slowly recognized that its essence was already analyzed by Friedrichs; 
references to an early work on the model can be found in |Exl notes to Sec. 3.2]. 
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The model exists in numerous modifications; we describe here the simplest one. 
We suppose that the state Hilbert space of the system has the form'H := C®iy^(M+) 
where the one-dimensional subspace is identified with the space mentioned 
above; the states are thus described by the pairs (^) with a G C and / G L^(M+). 
The Hamiltonian is the sclf-adjoint operator on H, or rather the family of self- 
adjoint operators labelled by the coupling constant 5 € M, defined by 

(3.1, H, = H„^,V. („t t/')' 

where Aq is a positive parameter, v e L^(K+) is sometimes called form factor, and 
Q is the operator of multiplication, {Qf){S,) = C.f (0- This in particular means that 
the continuous spectrum of Hq covers the positive real axis and the eigenvalue Aq is 
embedded in it; one expects that the perturbation gV can move the corresponding 
resolvent pole from the real axis to the complex plane. 

To see that it is indeed the case we have to find the reduced resolvent. The 
model is solvable in view of the Friedrichs condition, E^VE^ = where E^ is the 
projection to Hd L^(M^), which means that the continuum states do not interact 
mutually. Using the second resolvent formula and the commutativity of operators 
Etj and RHa{z) together with E^ + Ei = I we can write E^RHg{z)E-a as 

E^Rh,{z)E^ - gE^RH,{z)EuVE^RH^{z)E^ - gE^RH,{z)E^VEaRH,iz)E^ ; 

in a similar way we can express the 'off-diagonal' part of the resolvent as 

EaRh^{z)E^ = -gEdRH,{z)EdVE^RH^{z)E^ , 

where we have also employed the Friedrichs condition. Substituting from the last 
relation to the previous one and using {Hq — z)EuRho{z) = E^ together with the 
explicit form of the operators Hq, V , we find that R^^{z) acts for Im^ 7^ on 
■Hu = C as multiplication by the function 

(3.2) r^" : r^"(z) := (^-z + Aq + 5' dA^ . 

To make use of this result we need an assumption about the form factor, for instance 

(a) there is an entire / : C — >■ C such that [^^(A)^ = /(A) holds for all A G (0, 00) ; 

for the sake of notational simplicity one usually writes f{z) ~ \v{z)\^ for nonreal 
z too keeping in mind that it is a complex quantity. This allows us to construct 
analytic continuation of ?'g(-) over ac{Hg) = M+ to the lower complex halfplane in 
the form r{z) = \—z + w{z,g)\'^ , where 

w{\,g):=\,+g'l{\)-T^ig'\v{X)\^ ... A>0 

(3.3) 

w{z,g):= X^ + g'' r ^^^^Ai-2-Kig''\v{z)\^ ... Im^ < 
and /(A) is defined as the principal value of the integral, 

the analyticity can be checked using the cdge-of-thc-wedge theorem . 

These properties of the reduced resolvent make it possible to prove the mero- 
morphic structure of its analytic continuation. Resonances in the model are then 
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given by zeros of the function z i~> w{z,g) — z. An argument using the imphcit- 
function theorem |Ex[ sec. 3.2] leads to the following conclusion: 



Theorem 3.1. Assume (a) and v{Xo) 0, then r(-) has for all sufficiently 
small \g\ exactly one simple pole Zp{g) :— Xp{g) — i6p{g). The function Zp{-) is 
infinitely differentiable and the expansions 

(3.4) \p{g) = Xo+g'l{Xo)+0{g^), Sp{g) ^ 7rg'\v{Xo)f + 0{g^) , 

are valid in the vicinity of the point g = referring to the unperturbed Hamiltonian. 

To summarize the above reasoning we have seen that resonance poles can arise 
from perturbation of eigenvalues embedded in the continuous spectrum and that, 
at least locally, their distance from the real axis is the smaller the weaker is the 
perturbation. Moreover, one observes here the phenomenon called spectral con- 
centration: it is not difficult to check that the spectral projections of Hg to the 
intervals Ig (Aq — f3g, Xq — jig) with a fixed /3 > satisfy the relation 

s-liiji^ff,(/g) = E^. 

Friedrichs model also allows us to illustrate other typical features of resonant sys- 
tems. We have mentioned already the deep insight contained in the Fermi golden 
rule, which in the present context can be written as 

A=Ao 

where {E^^^} is the spectral decomposition of and Pc{Hq) the projection to the 
continuous spectral subspace of this operator. To realize that this is indeed what we 
known from quantum-mechanical textbooks, it is enough to realize that we use the 
convention h = 1 and formally it holds £^E^^^ Pc{Ho) = |A)(A|. Using the explicit 
form of the operators involved we find 

A — An 



dX Jo 



which is nothing else than the first nonzero term in the Taylor expansion (3.4). On 
the other hand, a formal use of the rule may turn its gold into brass: a warning 
example concerning the situation when the unperturbed eigenvalue is situated at 
the threshold of CTc(i?o) is due to J. Howland |Ho74j . see also |Ex[ Example 3.2.5]. 
Recent analysis of near-threshold effects in a generalized Friedrichs model together 
with a rich bibliography can be found in [D JNll] . 

Resonances discussed so far have been resolvent resonances. One can also 
consider the pair {Hg,Ho) as a scattering system. Existence and completeness of 
the wave operators is easy to establish since the perturbation gV has rank two. 
What we are interested in is the on-shell S-matrix: if v is piecewise continuous and 
bounded in IR+ one can check ^Ex., Prop. 3.2.6] that it acts as multiplication by 

S{X) = 1 + 2mg'^ lim_ \v{X)\^ r^{X + ie) . 

If V satisfies in addition the assumption (a) above, the function S{-) can be analyt- 
ically continued across R"*" . It is obvious that if such a continuation has a pole at a 
point z-p of the lower complex halfplane, the same is true for r(-), on the other hand, 
it may happen that a resolvent resonance is not a scattering resonance, namely if 
the continuation of has a zero at the point Zp. 
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Finally, the model can also describe a decaying system if we suppose that at 
the initial instant t — the state is described by the vector which span the one- 
dimensional subspace Hu- The main question here is whether one can approximate 
the reduced evolution by a semigroup in the sense of ( |2.4| ; a natural guess is that 
it should be possible in case of a weak coupling. Since the reduced resolvent is 



of the form (2.3) we can express the corresponding measure, calculate the reduced 



propagator according to ( 2.2 ) and estimate the influence of the second term in ( 2.4 ) 



This leads to the following conclusion, essentially due to |De76| : 

Theorem 3.2. Under the stated assumptions there are positive C, go such that 

\v{t) - Ae-''^'\ < ^ 

holds for allt>0 and \g\ < go with A := [1 — 5^/'(zp)]~^, where I{z) is the integral 
appearing in the second one of the formulcE \37[ 



The simple Friedrichs model described here has many extensions and in no way 
we intend to review and discuss them here limiting ourselves to a few brief remarks: 

(a) Some generalization of the model cast it into a more abstract setting, cf. for 
example |Mo96l IDR071 IDJNll] . Others are more 'realistic' regarding it as 
a description of a system interacting with a field, either a caricature one-mode 
one |DE87-89j or considerably closer to physical reality |BFS98L IHHH 08 in 
a sense returning the model to its Lee version which stimulated interest to it. 

(b) Friedrichs model clones typically use the simple procedure — attributed to 
Schur or Feshbach, and sometimes also to other people — we employed to 



get relation (3.2) expressing projection of the resolvent to the subspace 
sometimes it is combined with a complex scaling. 

(c) While most Friedrichs-type models concern perturbations of embedded eigen- 
values some go further. As an example, let us mention a caricature model of a 
crystal interacting with a field [DEH04j in which the unperturbed Hamilton- 
ian has a spectral band embedded in the continuous spectrum halfline referring 
to states of a lower band plus a field quantum. The perturbation turns the 
embedded band into a curve-shaped singularity in the lower complex halfplane 
with endpoints at the real axis. One can investigate in this framework decay 
of 'valence-band' states analogous to Theorem 3.2 etc. 

(d) The weak-coupling behavior described in Theorem 3.2 can be viewed also from 
a different point of view, namely that the decay law converges to a fixed expo- 
nential function as g when we pass to the rescaled time t' = g~^t. This is 
usually referred to as van Hove limit in recognition of the paper |vH55j : the 
first rigorous treatment of the limit belongs to E.B. Davics, cf. jPaj . 



4. Resonances from perturbed symmetry 

The previous section illustrates the most common mechanism of resonance 
emergence, namely perturbations of eigenvalues embedded in the continuum. A 
typical source of embedded eigenvalues is a symmetry of the system which prevents 
transitions from the corresponding localized state into a continuum one. Once such 
a symmetry is violated, resonances usually occur. Let us demonstrate that in a 
model describing a Schrodinger particle in a straight waveguide, perturbed by a 
potential or by a magnetic field, the idea of which belongs to J. Nockel |No92j . 
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4.1. Nockel model. We consider two-dimensional 'electrons' moving in a 
channel with a potential well. The guide is supposed to be either a hard-wall 
strip f2 := M X 5 with S = (—a, a), or alternatively the transverse confinement can 
be modelled by a potential in which case we have iS = M. The full Hamiltonian 
acting on H L^(f2) is given by 

(4.1) HiB, A) := (-za, - Byf + V{x) ~ + W{y) + \U{x, y) ; 

if is a strip of width 2a the transverse potential W may be absent and we 
impose Dirichlet conditions at the boundary, \y\ = a. The real- valued functions V 
describing the well in the waveguide — or a caricature quantum dot if you wish 
— and W are measurable, and the same is true for the potential perturbation 
U; further hypotheses will be given below. The number B is the intensity of the 
homogeneous magnetic field perpendicular to to which the system is exposed. 

The unperturbed Hamiltonian H{0) := H{0, 0) can be written in the form 
hX ® I + I ® which means its spectrum is the 'sum' of the corresponding 
component spectra. If the spectrum of the transverse part is discrete the 
embedded eigenvalues can naturally occur; we are going to see what happens with 
them under influence of the potential perturbation XU and/or the magnetic field. 
Let us first list the assumptions using the common notation (x) :— ^/l+x'^. 

(a) lim|2,|_j.oo W{x) — +oo holds if S* = M , 

(b) V ^0 and \V{x)\ < const (x)-'^-" for some e > 0, with V{x) dx < , 

(c) the potential V extends to a function analytic in Mao {z (z C : \ arg z| < ao} 
for some ao > and obeys there the bound of assumption (b), 

(d) \U{x,y)\ < const (x)^'^^^ holds for some e > and all {x,y) G H.. In addition, 
it does not factorize, U{x,y) ^ Ui{x)+U2{y), and U{-,y) extends for each fixed 
y € S to an analytic function in Aiao satisfying there the same bound. 

The assumption (a) ensures that the spectrum of := —dy + W{y), denoted 
as is discrete and simple, i^j+i > Vj. The same is true if = (— a, o) 

when we impose Dirichlet condition at y = ±a, naturally except the case when W 
grows fast enough as y — )• ±a to make the operator essentially self-adjoint. The 
assumption (b) says, in particular, that the local perturbation responsible for the 
occurrence of localized states is short-ranged and non-repulsive in the mean; it is 
well known that in this situation the longitudinal part hX :— + V{x) of the 
unperturbed Hamiltonian has a nonempty discrete spectrum, 

Hi < IJ.2 < ■ ■ ■ fJ-N < , 

which is simple and finite |Si76L IBGS77] ; the corresponding normalized eigenfunc- 
tions 0„, n = 1, . . . , iV , are exponentially decaying. 

To be able to treat the resonances we need to adopt analyticity hypotheses 
stated in assumptions (c) and (d). Note that in addition to the matter of our interest 
the system can also have 'intrinsic' resonances associated with the operator , 
however, the corresponding poles do not approach the real axis as the perturbation 
is switched off. In addition, they do not accumulate except possibly at the threshold 
}AC71j . and if V decays exponentially even that is excluded (Je78i. Lemma 3.4]. 

Since ac{h^) = [0, oo), the spectrum of the unperturbed Hamiltonian consists 
of the continuous part, fTc(^(0)) = crcss(-ff(0)) = [vi,oo), and the infinite family of 
eigenvalues 

ap(iJ(0)) = {fin + i^j ■■ n = l,...,N, j = l,2,...} . 
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A finite number of tliem are isoiated, wliile tlie remaining ones satisfying the con- 
dition /i„ + i>j > i^i are embedded in the continuum; let us suppose for simplicity 
that they coincide with none of the thresholds, /i„ + Vj ^ i>k for any k. 

To analyze the resonance problem it is useful to employ the transverse-mode 
decomposition, in other words, to replace the original PDE problem by a matrix 
ODE one. Using the transverse eigenfunctions, h^Xj = '^jXj > introduce the 
embeddings JT} and their adjoints acting as projections by 

J* : L^{n)^L\R), (J,*5)(x) = (x„g(x,.))i.(s) ; 
then we replace H{B, A) by the matrix differential operator {Hjk{B, A)}°°j,^j^ with 
H,k{B,X) := J*H{B,X)Jk ^ {^dl + V{x) + iy,)S,k+U,k{B,X), 



where mj^^ := Jg dy and Ujk{x) := JgU{x, y)Xj{y)Xk{y) dy. 

4.2. Resonances by complex scaling. Nockel model gives us an opportu- 
nity to illustrate how the complex scaling method mentioned in the introduction 
can be used in a concrete situation. We apply here the scaling transformation to 
the longitudinal variable starting from the unitary operator 

Sg: iSg^){x,y)^e'/^ij{e''x,y), 9eR, 

and extending this map analytically to Mao which is possible since the transformed 
Hamiltonians are of the form He{B, X) SeH{B, A)5e"^ = Hg{id) +Ue{B, A) with 

(4.2) He{0) -e-^'d^ ~ + Veix) + W{y) , 
where Vg{x) := V{e^^) and the interaction part 

Ug{B, A) := 2z e-'By 9, + S^y^ + \Ug{x, y) 

with Ug{x,y) :— U{e^^,y). Thus in view of the assumptions (c) and (d) they 
constitute a type (A) analytic family of m-sectorial operators in the sense of 'Ka] for 
|Im6'| < min{ao,7r/4} . Denoting Rg{z) := {Hg{0) - one can check ,DEMOl] 
that 

(4.3) \\Ug{B, X)Rg{iy, + Mi - 1)11 < c{\B\ + \B\^ + |A|) 

holds for jlm^l < min{ao, 7r/4}, and consequently, the operators Hg{B,X) also 
form a type (A) analytic family for B and A small enough. The free part (4.2) of 
the transformed operator separates variables, hence its spectrum is 

oo 

(4.4) a{Hg{0))^\J{i.,+a{h^)} , 

where :— —e^^^d^ + Vg{x). Since the potential is dilation-analytic by assump- 
tion, we have a typical picture: the essential spectrum is rotated into the lower 
halfplane revealing (fully or partly) the discrete spectrum of the non-selfadjoint 
operator hg which is independent of 6; we have 

(4.5) a (hn = e-2«M+ u {^i, • • ■ , Mjv} U {pi, p2 ■ • ■ } , 
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where pr are the 'intrinsic' resonances of hX . In view of the assumptions (c) and 
^■n + i^j 7^ for no k the supremum of Im over any finite region of the complex 
plane which does not contain any of the points Vk is negative, hence each eigenvalue 
+ hs-s a neighbourhood containing none of the points pk + Vj' ■ Consequently, 
the eigenvalues of Hg{0) become isolated once lm9 > 0. Using the relative bound- 



edness (4.31 we can draw a contour around an unperturbed eigenvalue and apply 
perturbation theory; for simplicity we shall consider only the non-degenerate case 
when fin + i^j =/= Pn' + i^j' for different pairs of indices. 

We fix an unperturbed eigenvalue cq — pn + Vj and choose 9 = ifi with a 



/? > 0; then in view of (4.4) and (4.5 1 we may chose a contour F in the resolvent 
set of He{0) which encircles just the eigenvalue cq. We use the symbol Pg for the 
eigenprojection of Hg{0) referring to eo and set 

^(,) 1 r Reiz) 



2'Ki Jp (eo — z)P 

for p = 0, 1, . . . , in particular, 5^°'' = —Pg and Re{z) Sg^'^ is the reduced resolvent 
of Hg{0) at the point z. The bound ( |4.3[ ) implies easily 

\\UgiB,X)Si''^\\ < (dist(r,eo))'^(|B| + \B\^ + |A|) 

with some constant c for all ImO G (0, ao) and p > 0. It allows us to write the 
perturbation expansion. Since eo = pn + fj holds by assumption for a unique pair 
of the indices, we obtain using fKaj Sec. II. 2] the following convergent series 



CO 



(4.6) e{B, A) = pn + Vj + e,„(B,A), 

m— 1 

where e^(i3, A) = Ep,+...+p„=„^-l Tr 0" i Ug{B, \)S'f^\ Using the above 

estimate we can estimate the order of each term with respect to the parameters. We 
find em{B,\) = EzHqC' (S'A"-'), in particular, we have e„(5) = 0(B'"), and 
em(A) = ©(A™) for pure magnetic and pure potential perturbations, respectively. 



The lowest-order terms in the expansion (4.61 can be computed explicitly. In 
the non-degenerate case, dimPg — 1, we have e{'^{B,X) ~ Tr {l4g{B, X)Pg). After 
a short calculation we can rewrite the expression at the right-hand side in the form 
'2-iBmfj' ((^„,0JJ + B^mfj' + A (0„, C/jj^„). Moreover, i (f>n) = {<i}n,idx4'n) is 
(up to a sign) the group velocity of the wavepacket, which is zero in a stationary 
state; recall that eigenfunction 0„ of is real-valued up to a phase factor. In 
other words, 

(4.7) er(i?. A) ^B^ [ \xMf dy + X [ U{x,y) \M^)xMf dxdy 
Js JRxS 

with the magnetic part independent of n. As usual in such situations the first-order 

correction is real-valued and thus does not contribute to the resonance width. 



The second term in the expansion (4.6) can be computed in the standard way 



|RS[ Sec. XII. 6]; taking the limit Im6' — > in the obtained expression we get 

oo 

(4.8) ef'(S,A) = -Y,{^MB,>^)(f>n, [{h'' -e„ + i,k-iOy'Ju,k{B, X)cf>n) . 

We shall calculate the imaginary part which determines the resonance width in the 
leading order. First we note that it can be in fact expressed as a finite sum. Indeed, 
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fceo := inax{fc : eo — I'k > 0} is finite and nonzero if tlie eigenvalue eo is embedded, 
otherwise we set it equal to zero. It is obvious that TZk ( {h^ — ep + J^fe~i0) ) 



is Hermitean for k > fce„, hence the corresponding terms in (4.8 1 are real and 

Ime^2"(B,A) = ^(W,fe(B,A)0„,(Im7efc)Z^,fe(S,A)(/.„) . 
fc=i 

The operators ImTZk can be expressed by a straightforward computation |DEMOl] . 
To write the resuh we need := [l+\V\^/'^{~dl-z)-'^\V\'^/'^sgn (V)] "\ in other 
words, the inverse to the operator acting as 

^ Jr 

We also need the trace operator : — > C which acts on the first Sobolev space 
VF^'^ as T^4> := (j){a^fE) for a — ± and E > Q where 4) is the Fourier transform 
of (j). Armed with these notions we can write the imaginary part of the resonance 
pole position up to higher-order terms as 

(4.9) ImefHS, A) = ^ E ^-7^= K-..w(eo-^^fc+iO) W,fe(i3, A)0„|' 

k 

fe=l<T = ± "^0 fc 

+2\B m!'^^ Im 



A2 

It. 

2 



(<o_^,w(eo-;^fe+iO)0'„ , <^_^^w(eo-J^fc+iO)?7jfe0„) 

,"„_,^w(eo-^/fe+^0) C/,fc0„|' } + 0(5^) + 0{B^\) , 
where as usual f{E + zO) = lim£^o+ f{E + ie). Let us summarize the results: 

Theorem 4.1. Assume (a)-(d) and suppose that an unperturbed eigenvalue 
eo = Mn + > ^1 simple and coincides with no threshold Vk ■ For small enough B 
and A the Nockel model Hamiltoni an J^.JP ha s a simple resonance pole the position 
of which is given by the relations |^.6[ )-(^.^). The leading order in the expansion 



obtained by neglecting the error terms in {4-9) is the Fermi golden rule in this case 



The symmetry in this example is somewhat hidden; it consists of the factorized 
form of the unperturbed Hamiltonian H(0) which makes it reducible by projections 
to subspaces associated with the transverse modes. It is obvious that both the 
potential perturbation — recall that we assumed U{x,y) ^ Ui{x)+U2{y) — and the 
magnetic field destroy this symmetry turning thus embedded eigenvalues coming 
from higher transverse modes into resonances. At the same time, the described 
decomposition may include other, more obvious symmetries. For instance, if the 
potential W is even with respect to the strip axis — including the case when 
S = {—a, a) and W = — the unperturbed Hamiltonian commutes with the 
transverse parity operator, 'ip(x,y) 1— )■ ip{x,—y), and the transversally odd states 
are orthogonal to the even ones so embedded eigenvalues arise. 

Nockel model is by far not the only example of this type. We limit ourselves here 
to quoting one more. Consider an acoustic waveguide in the form of a planar strip 
of width 2a into which we place an axially symmetric obstacle; the corresponding 
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Hamiltonian acts as Laplacian with Neumann condition at the boundary, both of 
the strip and the obstacle. Due to the axial symmetry the odd part of the operator 
gives rise to at least one eigenvalue in the interval (O, which is embedded 

into the continuous spectrum covering the whole positive real axis [ELV94| . If the 
obstacle is shifted by e in the direction perpendicular to the axis, such an eigenvalue 
turns again into a resonance for the position of which one can derive an expansion 
in powers of e analogous to Theorem 4.1 cf. [APVOO] . 



5. Point contacts 

The resonance models discussed in the previous two sections show that we 
should be more precise speaking about solvable models. The question naturally is 
what we have finally to solve when trying to get conclusions such as formulae for 
resonance pole positions. In both cases we have been able to derive weak-coupling 
expansions with explicit leading terms which could be regarded as confirmation of 
the Fermi golden rule for the particular model. One have to look, however, into 
which sort of problem the search for resonances was turned. For the Friedrichs 
model it was the functional equation w{z, g) = z with the left-hand side given by 



(3.3), and a similar claim is true for its clones, while in the Nockel model case we 
had to perform spectral analysis of the non-selfadjoint operatoi[^ 770(5, A). 

Not only the latter has been more difficult in the above discussion, the difference 
becomes even more apparent if we try to go beyond the weak-coupling approxima- 
tion. Following the pole trajectory over a large interval of coupling parameters 
may not be easy even if its position is determined by a functional equation and 
one have to resort usually to numerical methods, however, it is still much easier 
than to analyze a modification of the original spectral problem. Recall that for the 
Friedrichs model pole trajectories were investigated already in [Ho58j where it had 
been shown, in particular, that for strong enough coupling the pole may return to 
the (negative part of the) real axis becoming again a bound state. 

In the rest of this review we will deal with models which are 'solvable' at least 
in the sense of the Friedrichs model, that is, their resonances are found as roots 
of a functional — sometimes even algebraic — equation. In this section we will 
give examples showing that this is often the case in situations where the interaction 
responsible for occurrence of the resonances is of point or contact type. 



5.1. A simple two-channel model. The first model to consider here will 
describe a system the state space of which has two subspaces corresponding to two 
internal states; the coupling between them is of a contact nature. To be specific, 
one can think of a system consisting of a neutron and a nucleus having just two 
states, the ground state and an excited one. Their relative motion can be described 
in the Hilbert space L^(IR^) ® L^(IR'^); we suppose that the reduced masses is the 
two channels are the same and equal to ^ so that the Hamiltonian acts on functions 
supported away from the origin of the coordinates as — A and — A + £^, respectively, 
where i? > is the energy difference between the ground and the excited states. 

Before proceeding further, let us note that the above physical interpretation 
of the model coming from |Ex91| is not the only possible. The two channels can 



The same is true also for most 'realistic' descriptions of resonances using complex scaling, 
in particular, in the area of atomic and molecular physics — see, e.g., |RSI Sec. XII. 6] or ,Mo98] . 
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be alternatively associated, for instance, with two spin states; this version of the 
model was worked out in |CCF09] . also in dimensions one and two. 

To construct the Hamiltonian we start from the direct sum Aq = Aq.i ffi A0.2 
where the component operators act as Aq^i :— —A and Aq^2 ■— —A + E, respec- 
tively, being defined on W'^'^{M.'^ \ {0}). It is not difficult to check that ^0 is a 
symmetric operator with deficiency indices (2, 2); we will choose the model Hamil- 
tonian among its self-adjoint extensions. The analysis can be simplified using the 
rotational symmetry, since the components of Aq referring to nonzero values of the 
angular momentum are essentially self-adjoint, and therefore a nontrivial coupling 
is possible in the s-wave only. As usual wc pass to reduced radial wave functions 
/ : f{r) rtp{r) ; we take H = Hi ® H2 with Hj := L^(IR+) as the state space of 
the problem. The construction starts from the operator Hq = Hq i i?o.2, where 

ffoa:=-^, Ho,2:=-^+E, D{Ho,,) ^ W^^'iR+) , 

which has again deficiency indices (2, 2) and thus a four-parameter family of self- 
adjoint extensions. They can be characterized by means of boundary conditions: 
for each matrix -4 = ^) with a, 6 S M and c S C we denote by Hj\ the operator 
given by the same differential expression as Hq with the domain D{HyC) C D{Hq) — 
W^'^{R+) ® W'^''^{R+) specified by the conditions 

(5.1) /{(0) = a/i(0) + c/2(0), /^(0) = c/i(0) + 6/2(0); 

it is easy to check that any such is a self-adjoint extension of Hq. There may be 
other extensions, say, with decoupled channels corresponding to a = cx) or 6 = cxi 



but it is enough for us to consider 'most part' of them given by (5.1). 

If the matrix A is real the operator is invariant with respect to time re- 
versal. The channels are not coupled if c = ; in that case = Ha ® Hi, 
where the two operators correspond to the s-wave parts of the point-interaction 
Hamiltonians Ha.o and Hp^ in the two channels [AGHH] with the interaction 
strengths a := ^ and f3 :— respectively, and its spectrum is easily found. To 
determine cr(H_^) in the coupled case, we have to know its resolvent which can be 
determined by means of Krein's formula using the integral kernel GD{r,r';k) = 

diag( °''''''^'' 27fc'''''' '' ' ; 27k""''' '' ' )' where k := \/k^ — E, of the operator Hj^ 

with Dirichlet decoupled channels. The kernel of (Hj^ — z)^^ for z G p{HjC) equals 

GAr, r'; k) = G^{r, /; k) + D{k)-^ ^j^^^,,^ ^J^^^ j , 

where as usual k := and D{k) := {a — ik){b—iK) — |cp. 

It is straightforward to check that pole singularities of the above the resolvent 
can come only from zeros of the 'discriminant' D{k). In the decoupled case, i.e. if 
c = and Aq — (g , the expression factorizes, and consequently, it vanishes iff 
k — —ia or k = —ib. Several different situations may arise: 

• If a < the operator iJ^^ has eigenvalue —a^ corresponding to the eigenfunction 
f{r) = yj—1a ^ ) while for a > the pole now corresponds to a zero-energy 
resonance or an antibound state 



If 6 < 0, then Hj^^ has eigenvalue E — \? corresponding to /(r) = ^/—2b {^br), 



otherwise it has a zero-energy resonance or an antibound state. 
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The continuous spectrum of the decoupled operator covers the positive real axis 
being simple in [0,E) and of multiplicity two in [E,oo). We are interested mainly 
in the case when both a, b are negative and b^ < E ; under the last condition the 
eigenvalue of Hf, is embedded in the continuous spectrum of Ha- 

Let us next turn to the interacting case, c ^ 0. Since the deficiency indices of 
Hq are finite, the essential spectrum is not affected by the coupling. To find the 
eigenvalues and/or resonances of we have to solve the equation 



(5.2) (a-ik) [b-i^/k'^^E^ 



It reduces to a quartic equation, and can therefore be solved in terms of radicals; 
for simplicity we limit ourselves to the weak-coupling case when one can make the 
following conclusion |Ex91) . 

Theorem 5.1. (a) Let crp(i?^u) he simple, —a?' ^ E—b^, then the perturbed 
first-channel bound/ antihound state behaves for small \c\ as 

, , 9 2a|cP - E- bsj a^ + E , 

b + Va^+E Va2^{b+Va^+Ef 

In particular, zero- energy resonance corresponding to a — turns into an antibound 
state if Hj\^„ has an isolated eigenvalue in the second channel, b < —\fE, and into 
a bound state otherwise. 

(b) Under the same simplicity assumption, if Hj^^^ has isolated eigenvalue in the 
second channel, b < — \fE , the perturbation shifts it as follows 

a + Vb^-E ^/W^{a + VW^) 



On the other hand, if Hj\g has embedded eigenvalue, —\/E < b <0, it turns under 
the perturbation into a pole of the analytically continued resolvent with 

Ree.(c) = E-b^ + J^^+Oi\c\^), 



(c) Finally, let Hj\^„ have an isolated eigenvalue of multiplicity two, b — —\l o?^ E, 
then under the perturbation it splits into 

2a^+'ia^E+E^ 



ei,2(c) = -a^ T I4^a^fa^\c\ + ^^^^,^^^3^, |cp + 0{\c\^) . 

The model can be investigated also from the scattering point of view. Since 
the couplings is a rank- two perturbation of the free resolvent, the existence and 
completeness of the wave operators fi±(i?4, Uaq) follow from Birman-Kuroda the- 
orem (R,St Sec. XI. 3]. It is also easy to check that the scattering is asymptotically 
complete, what is more interesting is the explicit form of the S-matrix. To find it 
we look for generalized eigenfunctions of the form f(r) = (e"''^'' — Ae*'^'', Be'"'') 



which belong locally to the domain of Ha- Using boundary conditions (5.1 1 we find 



{a+ik){b-iK)-\c\^ 2ikc 

^ = ^°(') = W) ' W) 
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If the second channel is closed, < E, the scattering is elastic, \A\ ~ 1. We are 
interested primarily in the case when -ff^g has an embedded eigenvalue which turns 
under the perturbation into a resonant state whose lifetime is 

inspecting the phase shift we see that it has a jump by tt in the interval of width of 
the order 2 Im 62 (c) around Re 62 (c) . More specifically, writing the on-shell S-matrix 
conventionally through the phase shift as So{k) = e^*''"*^'') we have 

Oo(k) — arctan , (mod tt) . 

^ ' a{b + VE^) - \c\^ ^ ' 

The resonance is then seen as a local change of the transmission probability (and 
related quantities such as the scattering cross section), the sharper it is the closer 
the pole is to the real axis. This is probably the most common way in which 
resonances are manifested, employed in papers too numerous to be quoted here. 

On the other hand, if the second channel is open, fc^ > E, the reflection and 
transmission amplitudes given above satisfy + | \B\^ — 1. The elastic scattering 
is now non-unitary since B ^ which means that the 'nucleus' may now leave the 
interaction region in the excited state. The said relation between the amplitudes 
can alternatively be written as |5'o,i-).i(fc)P + |5'o,i^2(fc)P = 1 which is a part of 
the full two-channel S-matrix unitarity condition. 

The model also allows us to follow the time evolution of the resonant state. 



in particular, to analyze the pole approximation (2.4) in this particular case. The 
natural choice for the 'compound nucleus' wave function is the eigenstate of the 
unperturbed Hamiltonian, 

(5.3) / : fir) = V^b JJ,, 

Using the explicit expression of the resolvent we find 

(f (H k^)-^f)= \c\^ + {a~'ik)ib+iK) 
^•''^ ' ^' {b+iKY[\c\'^ ~ {a-ik){b-iK)\ 



The reduced evolution is given by (2.2); using the last formula, evaluating the 



integral by means of the residue theorem and estimating the remainder we arrive 
after a straightforward computation to the following conclusion |Ex91| . 

Theorem 5.2. Assume a 7^ and —\fE < b < 0. The reduced propagator of 



the resonant state {5.3) is given by 
46 



2(|a|-a)6 -^k^t , ^-^k^,t 



■ e 



(a'^^b^+E)^ ^/E^{a-iyfE-S^Y 



TT Jq {z^+ia^){z^-i{E-b^))^ 
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and the decay law is 



PA{t) 



,2(Im e2)t 



+ 



2|cMRe 



2{\a\-a)h 



ib 



1^ Uklt-T^/i) 



„2(Im 62)4 



z^e ^ * dz 



{z'^+ia'^){z^ -i{E-b^)f 



l + Odcp)) 



where e,- 



,{c) 



^3 



j = 1,2, are specified in Theorem 5.1 



Hence we have explicit formula for deviations from the exponential decay law. 
Some of its properties, however, may not be fully obvious. For instance, the initial 
decay rate vanishes, ^4(0+) = 0, since Im(/, {H^—X)~^f) = ©(A"^/^) as A — > 00, 
cf. Proposition |6.1| below. On the other hand, the long-time behavior depends 
substantially on the spectrum of the unperturbed Hamiltonian. If H_Ag has an 
eigenvalue in the first channel, then the decay law contains a term of order of |c|^, 
which does not vanish as t 00; it comes from the component of the first-channel 



bound state contained in the resonant state (5.3) 



The last name fact it is useful to keep in mind when we speak about an unstable 
state lifetime. It is a common habit, motiveted by the approximation (2.4 1, to 



identify the latter with the inverse distance of the pole from the real axis as we 
did above when writing T'(c). If the decay law is differentiable, however, —P^{t) 
expresses the probability density of decay at the instant t and a simple integration 
by parts allows us to express the average time for which the initial state survives 
as r4 = PA{t) di; this quantity naturally diverges if limf_j.oo PA{t) 7^ 0. 

5.2. K-shell capture model: comparison to stochastic mechanics. The 

above model has many modifications, we will describe briefly two of them. The first 
describes a /3-decay process in which an atomic electron is absorbed by the nucleus 
and decays through the reaction e+p n+v with a neutrino emitted. One usually 
speaks about a K-shell capture because the electron comes most often from the 
lowest energy orbital, however, from the theoretical point any orbital mode can be 
considered. We assume again spherical symmetry and take % = %i ®%2 with 
Hj := L^(M+) as the state space. The first component refers to the (s-wave part 
of) atomic wave function, the other is a caricature description of the decayed states; 
we neglect the fact that neutrino is a relativistic particle. 

The departing point of the construction is again a non-self adjoint operator of 
the form Hq = Ho i ffl Ho_2, the components of which act as 



(5.4) 



1 d2 



2m dr^ 

with the domains D(i?o,i) = {f & W^'^ 



Vi{r) , Ho,2 ■■= - 



1 d2 



- E 



2M dx^ 

u{Rj-) = u{Rj+) = 0} for fixed 
< i?i < i?2 < • • ■ < Rn — we add the requirement /(0+) = if the angular 
momentum l = Q — and D{Ho,2) = {W^'^{R+) : /(0+) = /'(0+) = 0}. Here as 
usual VeXr) = V{r) + ^2mr2 and the potential is supposed to satisfy the conditions 



lim V{r) 



0, 



lim sup rV{r) 



7 < 00 , 
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under which the operator Hq is symmetric with deficiency indices (N + 1, iV + 1). 
Of all its self-adjoint extensions we choose a subclass that (i) allows us to switch 
off the coupling and (ii) couples each sphere locally to the other space. The adjoint 
operator Hq acts again as (5.4); the extensions H{a) described by a = (ai, . . . , a^) 
are specified by the boundary conditions 

N 

u'{R,+)-u'{Rj-)^ajf{0+), j^l,...,N and ^ a,u(i?j) - - ^ /(0+) ; 

i=i 

it is easy to see that under them the appropriate boundary form vanishes, the chan- 
nels are decoupled for a — 0, and the Hamiltonians Ha are time-reversal invariant. 

To solve the resonance problem in the model we have to find the resolvent of 
H{a) which can be again done using Krein's formula. We will describe the resolvent 
projection to the 'atomic' channel referring for the full expression and the proof to 
|ET92j . We introduce the kernel 

Gi(r,s;fc2) = — — ^Wfc(?'<)wfe(r>), 

W(vk,Uk) 

where as usual r< :— min{r, s}, r> :— max{r, s}, and the functions Uk,Vk are 
solutions to u" + Vu — such that Ufc(O-l-) = and is L'^ around cx), 

and furthermore, VF(w/c,Wfe) := ffc(r)u'j,(r) — v'f,{r)uk{r) is their Wronskian. 

Before stating the result, let us mention that the model can also cover the 
situation when the electron can be absorbed anywhere within the volume of the 
nucleus approximating this behavior by a family of equidistant spheres with Rj :— 
jR/N, j = 1,...,N, where R is the nucleon radius. Let a : [0,i?] — > M be a 
bounded piecewise continuous function understood as 'decay density', and take 
i?(a^) corresponding to aj^ :— ^a('^). On a formal level, the limit N ^ oo 
leads to an operator describing the two channels coupled through the boundary 
conditions 

u"{r) = a(r)/'(0+) , a{s)u{s) ds = - /(0+) , 

however, we limit ourselves to checking the strong resolvent convergence |ET92j . 

Theorem 5.3. The projection of the resolvent {H{a) — z)^^ to the 'atomic' 
channel is an integral operator with the kernel 

>^ iKMajakGi{r, Rj-, z)Gi{Rk, s; z) 
Gi[r,s;z)+ ^ — — tv ~T7rr^~^T7\- 



iKMJ2i 1=1 aiaiGi{Ri, Ri] z 



The projections of [H{a^) — z) ^ converge as N —i' oo to operator with the kernel 

G,{r,s-z)+ i^-M<t>.{r)Ms) 

m — iKil/ Jq a(r)a(s)Gi(r, s; z) drds 

where <j)k '■— a{s)Gi{-, s; z) ds, in the strong resolvent sense. 

The singularities correspond to zeros of the denominators in the above expres- 
sion. As an example, consider the 'atom' with Coulomb potential. 
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which has in the decoupled case, a = 0, poles at fc„ — n — 1,2,.... The 

Green function Gi can be expressed in terms of the standard Coulomb wave func- 
tions ipnimir,!!}, if) = Rni{r)Yira{'& , ■ To analyze weak-coupling behavior of the 
poles, we restrict ourselves to two cases, a surface-supported decay when = 1 
and a volume-supported decay with a constant 'density' when a(r) = a is a constant 

function on [0, i?]. We put k„ :— \J 2mE — (^)^ and introduce the form factor 



Bn{r) 



RRni{r) ... surface-supported decay 

r Rne{r) dr ... volume-supported decay 



A straightforward calculation |ET92j then yields the shifted pole positions, 
2m 2n 4 

We are interested particularly in the situation where the unperturbed eigenvalue is 

embedded, n > ^ ^^^j when k„ is real and the coupling shifts the pole into the 
lower complex halfplane giving rise to the resonant state with the lifetime 

In the case of the real decay, of course, all the unperturbed eigenvalues are embedded 
and the K-shell contribution is dominating. It has the shortest lifetime since i?„£(0) 
is nonzero for £ = only and m-jR <C 1, typically of order 10^^, so the form factor 
value is essentially determined by the wave function value at the origin. 

The K-shell capture model allows us to make an important reflection concern- 
ing relations between quantum and stochastic mechanics. The two theories are 
sometimes claim to lead to the same results [Ne] and there are cases when such a 
claim can be verified. The present model shows that in general there is a principal 
difference between the two. One can model such a decay in stochastic mechanics 
too considering random electron trajectories and summing the decay probabilities 
for their parts situated within the nucleus. The formula is given in |ET92I and 
we are not going to reproduce it here; what is important that the total probability 
is expressed as the sum of probabilities of all the contributing processes. In the 
quantum-mechanical model discussed here, on the other hand, one adds the ampli- 
tudes — it is obvious from the form factor expression in case of a volume-supported 
decay — and the total probability is the squared modulus of the sum. 

5.3. A model of heavy quarkonia decay. Let us finally mention one more 
modification of the model, this time aiming at description of decays of charmonium 
or bottomium, which are bound states of heavy quark-antiquark pairs, into a meson- 
antimeson pair. Such processes are known to be essentially non-relativistic; as an 
example one can take the decay ■0"(377O) — ^ DD where the D meson mass is 
« 1865MeV/c , thus rest energy of the meson pair is two orders of magnitude 
larger than the kinetic one released in the decay. 

If the interaction responsible for the decay is switched off the quark and meson 
pairs are described by the operators 

Hqj := - ^Aji - ^Aj2 -I- Vj{\xji-Xj2\) + 2mjc^ , 
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where mi and m2 are the quark and meson masses, respectively. As before we 
separate the center-of-mass motion and use the rotational invariance. Adjusting 
the energy threshold to 2m2C^ we can reduce the problem to investigation of self- 
adjoint extensions of the operator TJq^-* := H^^l © Hq^I on L^{Ri, oo) L^(R2, oo) 
defined by 

^"■^ - " ^ dTf + ^^"("^'^ + + ^^"^^ 

with D^Hq^j) :— C^(i?j, oo). Let us hst the assumptions. We suppose that the 
quarks can annihilate only they 'hit each other', i?i = 0, while for mesons we allow 
existence of a hard core, R2 = R > 0. One the other hand, the mesons are supposed 
to be non-interacting, V2 = 0; this may not be realistic if they are charged but it 
simplifies the treatment. In contrast to that, the interquark potential is confining, 
limr-j-oo Viir) = 00; we also assume that Vi € Lj^^ and a finite limr_>.o+ Vi{r) exists. 

To couple the two channels the deficiency indices of Hq^^ have to be (2,2); 
since we have put i?i = it happens only if £ = and we drop thus the index £ 
in the following. We will not strive again to describe all the extensions and choose 
a particular one-parameter family: the domain of the extension Ha will consist of 
functions / G VK^'^(M"'") ® W^'^{R, 00) satisfying the conditions 

(5.5) /i(0) = a/^(i?), /2(i?) = !!^a/{(0), 

mi 

with a G C. In the decoupled case, a = 0, we get Dirichlet boundary condition in 
both channels as expected; for a G M the Hamiltonian is time-reversal invariant. 

As the first thing we have to find the resolvent of iJ^, in particular its projection 
to the quark channel. In analogy with the previous section we can write its integral 
kernel Go{r, s; z) in terms of two solutions of the equation 

( - — + Vi{r) + 2(mi - m2)c2)/(r) = z/(r) 
V raj ar'^ / 

for z ^ M such that u(0) = and v is at infinity. Krein's formula helps again; 
by a straightforward computation AESS94] we get 

-ikm2\a\'^v{r]z)v{s]z) 

Ga{r, s; z) = Go r, s; z) H — , 

miv[0; z)U(v, a; z) 

where the denominator is given by 

D{v, a; z) := t;(0; z) — ifclaP — ^ w'(0; z) . 

mi 

The singularities are again determined by zeros of the last expression. One can work 
out examples such a natural confining potential, Vi{r) = ar+VQ+2{m2—mi)c^ , and 
its modifications, in which the resonance width can be expressed through the value 
of the quark wave function at the origin. This appears to be the case generally. 



Theorem 5.4. |AESS94j Under the stated assumptions, the quarkonium decay 
width is given for the n-th s-wave state by 

(5.6) r„(a) =87rfc„^|ap|V'„(0)p+O(|aH, k^ VmiE„ , 

mj 

provided the bound-state energy E^, adjusted by the difference of the rest energies, 
is positive; V'n(O) is the value of the corresponding wave function at the origin. 
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Note that while we have assumed the quark potential to be below bounded at 
the origin, the assumption can be relaxed. The theorem holds also for potentials 
with sufficiently weak singularity, in particular, for the physically interesting case 
of a linear confinement combined with a Coulomb potential. 



6. More about the decay laws 

Let us return to the time evolution of unstable systems, in particular, to proper- 
ties of the decay laws. In addition to the elementary properties mentioned together 



with the definition (2.11 we know so far only that in the weak-coupling situation 
they do not differ much from an exponential function coming from the leading term 
of the pole approximation. This says nothing about local properties of the decay 
laws which is the topic we are going to investigate in this section. 

Historically the first consequence of non-exponentiality associated with the be- 
low bounded energy spectrum concerned the long-time behavior of the decay laws; 
already in Kh57^ it was observed that a sharp energy cut-off leads to the C'(t~'^/^) 
behavior as t — ?► oo, and other examples of that type followed. Moreover, it is even 
possible that a part of the initial state survives the decay; we have seen a simple 



example at the end of Sec. 5.1 and another one will be given in Sec. [9]below. Here 



we concentrate on two other local properties of decay laws. 



6.1. Initial decay rate and its implications. The first one concerns the 
behavior of the system immediately after its preparation. Exponential decay has a 
constant decay rate which, in particular, means it is nonzero at t = 0. This may 
not be true for other decay laws. We note, for example, that is by definition 
an even function of t, hence if the (two-sided) derivative Pjp{0) exists it has to be 
zero. This happens for vectors from the form domain of the Hamiltonian: we have 
KV', e~*^*V')P < Pip{t) ^ 1 which leads easily to the following conclusion [HE73j . 

Proposition 6.1. Ifipe Q{H) the decay law satisfies P^(0+) = 0. 

The importance of this result stems from the peculiar behavior of unstable sys- 
tems subject to frequently repeated measurements knows as quantum Zeno effect, 
namely that in the limit of permanent measurement the system cannot decay. This 
fact was known essentially already to von Neumann and Turing, in the context 
of unstable particle decay it was first described by Beskow and Nilsson |BN67| 
followed by a serious mathematical work |Fr72|, ICh| which elucidated the mech- 
anism. It became truly popular, however, only after the flashy name referring to 
Zeno's aporia about a flying arrow was proposed in |MS77| . Since then the effect 
was a subject of numerous investigations, in part because it became interesting 
also from experimental and application points of view. However, since Zeno-type 
problems are not the subject of this survey we limit ourselves to quoting the review 
papers j Sch04|, IFP08| as a guide to further reading, and will discuss the topic only 
inasmuch it concerns the initial decay rate. 

Suppose that we perform on an unstable system a series of measurements at 
times t/n, 2t/n . . . , in which we ascertain that it is still undecayed. If the out- 
come of each of them is positive, the state reduction returns the state vector into 
the subspace "Hu and the resulting non-decay probability is 



M„(t) = P^{t/n)P^^{t/n) ■ ■■P^^_^{t/n) , 
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where ipj+i is the normalized projection of e on T-L^ and ijjQ := ijj, in 

particular, Mn{t) — {P^,{t/n))" if diniHu = 1 and all the vectors tpj coincide with 
up to a phase factor. Since lim„_j.oo(/(i/'^)" = exp{— /(0+)t} holds whenever 
/(O) = 1 and the one-sided derivative /(0-I-) exists, we see that P^(O-I-) = implies 
the Zeno effect, M{t) :— lim„^oo Mn{t) = 1 for all i > 0, and the same is true if 
dim^u > 1 provided the derivative P^(O-I-) has such a property for any S Hu- 

At the same time the above simple argument suggests that an opposite situa- 
tion, an anti-Zeno effect, is possible when P^(0-f ) is negative infinite; then M{t) = 
for any i > which means that the decay is accelerated and the unstable system 
disappears once the measurement started. The possibility of such a behavior was 
mentioned early |CSM77] . however, the attention to it is of a recent date only — 
we refer again to the review work quoted above. Before proceeding further we have 
to say that the two effects are understood differently in different communities. For 
experimental physicists the important question is the change of the observed life- 
time when the measurement are performed with a certain frequency, on the other 
hand a theoretical or mathematical physicist typically asks what happens if the 
period between two successive measurements tends to zero. 

Let us return to the initial decay rate. It is clear we have to estimate 1— for 
small values of t, which we can write as 2 Re {^|J, Eu{I—e~^^'^)^p) — \\Ea{I—e~^^*)^\\^, 
or alternatively cast it using spectral theorem into the form 



4 |°°sin2^d||i?fV^|p-4 |_ 



where the non-decreasing projection-valued function A i— )• := Eh{{—oo,X]) 
generates spectral measure Eh of the Hamiltonian H. By Schwarz inequality the 
above expression is non-negative; we want to find tighter upper and lower bounds. 

To this aim we choose an orthonormal basis {xj} in the unstable system sub- 
space Hu and expand the initial state vector as -0 = J2j '^jXj with |cjp = 1. 
The second term in the above expression can be then written as 



-4E E 



xt 



^^*/^ sin^d(x™,i?fx.) 



where dwjfc(A) := A.{xj,E^Xk) are real- valued measures symmetric with respect to 
interchange of the indices. Since the measure appearing in the first term can be 
written as d||i?|'^V''lP — J2jk CjCkdi^jk{X), the decay probability becomes 



sin —dwjk{X) 



l-PAt)=4j2c,CklJ_ 

jk 

-E /j"^'*^' sin^d..,,„(A) l_ 



if dim?^u = oo the involved series can easily be seen to converge using Parseval 
relation. Using next the normalization dujjkiX) — Sjk we arrive after a simple 
calculation |Ex05| at the formula 



(6.1) 



/OO n 
/ 
jKrn -ooJ-oo 



gin2 J^^* dujjrniX)dLOkm{p-) 



22 



PAVEL EXNER 



Consider first an upper bound. We fix a e (0,2] and use the inequalities \x\" > 
I sinx|" > sin^ x together with |A - < 2"(|A|" + |^|") to infer that 

1 — P/ (f'] f°° f°° 

VI ) < 2i-" V c.cfe / / |A - Ml" dio,^{X)dujkM 



J km 



OO 



/OO poo 
/ (iAr + i/^r) dw,™(A)dwfc™(M) < 4(|i/r)^ , 
jKm -ooJ-oo 

which means that 1 - P^{t) = 0(P) if i/) e Dom If this is true for some 

a > 1 we get P^{0+) = which a shghtly weaker result than Proposition [o] Note 
also that if dim — 1 and the spectrum of H is absolutely continuous there is an 
alternative way to justify the claim using Lipschitz regularity, since P{t) ~ \uj(t)\'^ in 
this case and Cu is bounded and uniformly a-Lipschitz iff a;(A) (1 + | A|") dA < oo. 

A lower bound is more subtle. We use the inequality I sin | > C\X — /i|t 



which holds with a suitable C > for A, G [— l/t, to estimate (6.1) as follows 

»i/t /.l/t 



1 - P^ it) > 2Ch^ I I (A - m)' (V-, dE^E^dE^i,) 
J -l/t J -l/t 



where Hn denotes the cut-off Hamiltonian, HEh{^n) with Ajv {—N,N). Di- 
viding the expression at the right-hand side by t and choosing t — N^^, we arrive 
at the following conclusion. 

Proposition 6.2. The initial decay rate of e Hu satisfies P^{0+) = — oo 
provided {{H%E^EHiAN))i, - \\PHNip\\^y^ = o{N) holds as N ^ oo. 

To illustrate how does the initial decay rate depend on spectral properties of 
the decaying state, consider an example in which dimT^u = 1, the Hamiltonian is 
bounded from below and -0 from its absolutely continuous spectral subspace is such 
that d(i/', E^ip) = uj{X) dX wher e w(A) « cX~'^ as A — >■ -|-oo for some c > and 



/3 > 1. If /3 > 2, Proposition 6.1 implies P^(O-t-) = 0. On the other hand, one can 



J^j^uj(X) dX tends to one, while the other two integrals diverge giving 



easily find the asymptotic behavior of the quantity appearing in Proposition 16.21 

N 

-JV' 



N rN / i-N /„\2 



-N 




A2 d^(A) / dL^(A) - / Ad^(A) « ^— N^-^ - — - N^-^f^ 



3-/3 



as +00, and consequently, P{Q+) = — oo holds for /3 e (1, 2). This shows that 



the exponential decay — which requires, of course, (j{H) = M by Theorem 2.1 — 
walks a thin line between the two extreme initial-decay-rate possibilities. Let us 
remark finally that while 'Zeno' limit is trivial for the exponential decay, it may 
not exist in other cases with /3 = 2; in [Ex, Rem. 2.4.9] the reader can find an 
example of such a distribution with a sharp cut-off leading to rapid oscillations of 
the function t H> (V',e~*^*V') which obscure existence of the limit. 
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6.2. Irregular decay: example of Winter model. Now we turn to another 
decay law property. In the literature it is usually tacitly assumed that P^{-) is a 
'nice', i.e. sufficiently regular function, typically by dealing with its derivatives. 
Our aim is to show that this property cannot be taken for granted which we are 
going to illustrate on another well-known solvable model of decay. 

An inspiration comes from the striking behavior of some wave functions in a 
one-dimensional hard- wall potential well observed in |Be96l ITh) . The simplest 
example concerns the situation when the initial function is constant (and thus 
not belonging to the domain of the Dirichlet Laplacian): it evolves into a steplike 
4'{x, t) for times which are rational multiples of the period, t ~ qT with q — N/M, 
and the number of steps increases with growing M, while for an irrational q the 
function 'ijj{x, t) is fractal with respect to the variable x. One may expect that such 
a behavior will not disappear completely if the hard wall is replaced by a singular 
potential barrier. It was illustrated in a 'double well' system |VDS02] : here we 
instead let the initial state decay into continuum through the tunneling. 

Decay due to a barrier tunneling is among the core problem of quantum me- 
chanics which can be traced back to Gamow's paper |Ga28| . The model in which 
the barrier is a spherical (5-shell is usually referred to as Winter model after the 
paper |Wi61| where it was introduced. A thorough analysis of this model can be 
found in |AGS87] : it has also various generalizations, we refer to [AGHH] for a 
bibliography. The Hamiltonian acting in is of the form 

Ha = -A + a5{\r\- R), a>0, 

with a fixed i? > 0; as usual we employ rational units, h — 2m — 1. For simplicity 
we restrict our attention to the s-wave part of the problem, using the reduced wave 
functions ^/j{f,t) = -^^r~'^(f>{r, t) and the corresponding Hamiltonian part, 

d2 

ha — -"T^ + aSir — R) ] 
dr-' 

we are interested in the time evolution, ip{f,t) = e~*^°*'(/'(r, 0) for a fixed initial 
condition ■(/;(r, 0) with the support inside the ball of radius R, and the corresponding 

decay law P^{t) = J^^ \(t>{r, t)\'^ dr referring to Hu = 1.^(5^(0)). 

It is straightforward to check |AGS87] that Ha has no bound states, on the 
other hand, it has infinitely many resonances with the widths increasing logarith- 
mically with respect to the resonance index [EF06J; a natural idea is to employ 
them as a tool to expand the quantities of interest |GMM9"5] . In order to express 
reduced evolution in the way described in Sec. [2] we need to know Green's function 
of the Hamiltonian ha which can be obtained from Krein's formula, 

{ha - e)-\r, r') = {ho - k^)-\r, r') + A(fc)a>,(r)$fe(/) , 

where ^^(r) := GQ{r, R) is the free Green function with one argument fixed, in 
particular, ^^(r) = ^ sin(fcr) e*''^ holds for r < R, and A(fc) is determined by 
J-interaction matching conditions at r = i?; by a direct calculation one finds 
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Using it we can write the integral kernel of e as Fourier transformation, 
u{t, r, r') = p{k, r, r')e~*'^ * 2k dk, where the explicit form of the resolvent gives 

, 2fcsin(fcr) sin(/cr') 

'^'^ ' ~ 7r(2fc2 + 2a2 sin^ kR + 2ka sin 2kR) ' 

The resonances understood as poles of the resolvent continued to the lower halfplane 
appear in pairs, those in the fourth quadrant, denoted as kn in the increasing order 
of their real parts, and —kn] we denote S — {kn, —kn, kn, —kn : n g N}. In the 
vicinity of fc„ the function p{-, r, r') can be written as 

II. i\ * Vn{r)vn{r') 

where w„(r) solves the differential equation haVn{r) = kf^Vni^) and x is locally ana- 
lytic. It is not difficult to see that the function p{-, r, r') decreases in every direction 
of the fc-plane, hence it can be expressed as the sum over the pole singularities, 

1 

k-k 

kes 

and the residue theorem implies X^fces ReSj,p(fc, r, r') = 0. Using these relations 
and denoting fc_„ := —kn with u_„ being the associated solution of the equation 
HaV^n{f) — k'^n'^-nif), wc arrivc after a short computation |EF07| at 



p{k, r, r') = ^ ^ ReSj, p{k, r, r') 



u{t,r,r') = M{kn,t)vn{r)vn{r') 



with M{kn,t) — \ e"" erfc(u„) and u„ := — e ^'^l'^kn\/t, leading to the decay law 

P^^{t) = CnClInlM{kn,t)M(kut) 

with Cn ■= 4'{r,0)vn{r) dr and /„; :— J^Vn{r)vi{r)dr; in our particular case 
we have Vn{r) — \f2Qn sin(fc„r) with the coefficient Qn equal to 

-2^1 ^ 



2kn + a^Rsui2knR + ashi2knR + 2fc„Q!i?cos 2fc„_R^ 

These explicit formula allow us to find (t) numerically. Let us quote an example 
worked out in |EF07j in which R — 1 and a = 500; the initial wave function is 
chosen to be constant, i.e. the ground state of the Neumann Laplacian in L^(B]^{0)) 
which corresponds to 4'{r,0) = R^^^^ \/^r X[o,h]{''') ■ 

The respective decay law is plotted in Figure [T] we see that it is irregular 
having 'steps', the most pronounced at the period T = 2R? /t: and its simple rational 
multiples. This is made even more visible from the plot of its logarithmic derivative 
(for numerical reasons it is locally smeared, otherwise the picture would be a fuzzy 
band). It is reasonable to conjecture that the function is in fact fractal. 

Let us add a few heuristic considerations in favor of this conjecture concerning 
the behavior of the derivative in the limit a — ^ oo. We can write the wave function 
as 0(r, t) c„ exp(— zfc2t)w„(r) where resonance position expands for a fixed n 

around kn^ := rnr/R as fc„ « fc„^o — kn,oi'^R)~^ + knfl{aR)~'^ — ik"^ i^{a^R)^^ . In 

the leading order we have Vn{r) ~ \f\ sin(fc„r) and the substantial contribution to 

the expansion of (j){r,t) comes from terms with n < [a^^'^^] for some < e < 1/3. 
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0.94 I 1 > 1 1 

T/3 T/2 T 

t 

Figure 1. Decay law for the initial state 0) = i?^^/^-\/3r; 
the inset shows its logarithmic derivative averaged over intervals 
of the length approximately T/200. 



The derivative of the decay law P,p.a{-) can be identified with the probability 
current through the sphere, P^^a{t) = — 2Im i)). To use it we have to 

know the expressions on the right-hand side; using the above expansion we find 

^(^' ^) - ./|f:(-l)"c„e--^o*(-^)e-^* f_kr^_kl^\ 

n— 1 ^ ^ 

and a similar expansion for (j)'{R,t) with the last bracket replaced by fc„_o- We 

observe that ^^=1 exp (-25^*) « f iif^'^^' a^y^i^ holds for j > -1 

where on the right-hand side we have denoted Ij := e~^^ dx — |r (^^)- 

Using this result we can argue that the decay law regularity depends on the 
asymptotic behavior of the coefficients c„. Suppose for simplicity that it is power- 
like, c„ — ©(fc^p) as rt — > oo. If the decay is fast enough, p > 1, we find that 
I^V.qIOI — const Q!'*/'^^'*/'^P — ^ holds as a oo uniformly in the time variable. 
The situation is different if the decay is slow, p < 1. Consider the example men- 
tioned above leading to the decay law featured in Figure [l] where c„ — (— 1)"^^ 
Since the real parts of the resonance poles change with a, it is natural to look at 
the limit of P^^a{ta) as a — >■ cx) at the moving time value ta '■= t(l + 2/aR). 

For irrational multiples of T we use the fact |BG88| that the modulus of 
J2n=i '^^^^ * fo'^ irrational t bound by C L^^^ where C,e depend on t only. In 
combination with the estimate, J2'^=i o-nbn < J2'^=i I X]j=i \ 1^" ^ ^n+i| we find 

^e-^'=".o*e-W*/fc^^„ < const a2/30-+i-^), 
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and consequently, P^^a{ia) — > as a — > oo. Assume next that t — with 
p,q £N. If pq is odd then SL{t) X]n=i repeatedly retraces according to 

|BG88| the same pattern, hence P^^aita) — an illustration can be seen in 
Fig. [l] at the half period. On the other hand, for pq even |S'l(OI grows linearly with 
L, and consequently, liuia^ao Pip,aita) > 0. For instance, a direct computation 
|EF07] yields the value at the period, limQ^^oo P^,a{Ta) = ~ -0.77. 

As the last remark in this section, we note that there is a relation between a lack 
of local regularity of the decay law and the 'anti-Zeno' property of Proposition |6]2j 
both occur if the energy distribution of the decaying state has a slow enough decay 
at high energies. The connection is no doubt worth of further exploration. 

7. Quantum graphs 

Many quantum systems, both spontaneously emerging in Nature and resulting 
from an experimentalist's design, no doubt intelligent one, have complicated geo- 
metrical and topological structure which can be conveniently modeled as a graph 
to which the particle motion is confined. Such a concept was first developed for 
the purpose of quantum chemistry |RS53j . however, it became a subject of in- 
tense investigation only at the end of the 1980's when tiny graph-like structures 
of semiconductor and other materials gained a prominent position in experimental 
physics. The literature on quantum graphs is vast at present; we limit ourselves 
with referring to the proceedings volume [EKKSTj as a guide for further reading. 

Quantum graphs are usually rich in resonances; the reason, as we see below, is 
that their spectra often exhibit embedded eigenvalues which, as we know, are typ- 
ically sources of resonance effects. Before we turn to the review let us briefly men- 
tion that while describing real-world quantum system through graphs is certainly 
an idealization, they can be approximated by more realistic 'fat-graph' structures 
in a well-defined mathematical sense; from our point of view here it is important 
than such approximations also include convergence of resonances |EP07| . 

7.1. Basic notions. As a preliminary, let us recall some basic notions about 
quantum graph models we shall need in the following. For the purpose of this 
review, a graph F consists of a set of vertices V = {Xj : j € /}, a set of finite edge^ 
C — {Cjn ■ {Xj,Xn) £ Ic C / X /}, and a set of infinite edges, sometimes also called 
leads, Coo = {'^joo '■ Xj G /c} attached to them. We consider metric graphs which 
means that each edge of F is isomorphic to a line segment; the notions of finiteness 
or (semi)infiniteness refer to the length of those segments. As indicated, we regard 
F as the configuration space of a quantum system with the Hilbert space 

H = L2([0,?,])e i^([0,oo)), 

the elements of which are columns ^ ~ {{fj : Cj G £}, {gj : £joo G £00})"^- 

For most part of this section we will suppose that the motion on the graph edges 

.2 

is free, i.e. governed by the Hamiltonian which acts there as — with respect to 
the arc-length variable parametrizing the particular edge. In order to make it a 
self-adjoint operator, we have to impose appropriate boundary conditions which 

^We assume here implicitly that any two vertices are connected by not more than a single 
edge and that the graph has no loops, which is possible to do without loss of generality since we 
are always able to insert 'dummy' vertices into 'superfluous' edges. 
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Figure 2. The model Tq for a quantum graph T with N internal 
finite edges and M external leads 



couple the wave functions at the graph vertices. One of the possible general forms 
of such conditions jGGl IHaOOL IKSOOj is 

(7.1) iU,-I)^,+ziU,+I)^^0, 

where Uj are unitary matrices, and and are vectors of the functional values 
and of the (outward) derivatives at the vertex Xj; in other words, the domain of 
the Hamiltonian consists of all functions on T which are locally W'^''^ and satisfy 
conditions (7.1 ). Note that coupling is local connecting boundary values in Xj only. 



Since handling Hamiltonians of graphs with a complicated topology may be 
cumbersome, one can employ a trick proposed in [Ku08j replacing T with the graph 
Tq in which all edge ends meet in a single vertex as sketched in Fig. |2] the actual 
topology of r will be then encoded into the matrix which describes the coupling in 
the vertex. Denoting TV = jj £ and M = Crx, we introduce the {2N + M)-dimen- 
sional vector of functional values by VE* = (^'f , . . . , ^«v)^ ^^'^ similarly the vector 



of derivatives \E'' at the vertex; the conditions (7.1) can be concisely rewritten as 
coupling on Tq characterized by (2iV + M) x (2N + M) unitary block-diagonal 
matrix U, consisting of the blocks Uj, in the form 

(7.2) (f/ - 7)^' + i(f7 + /)*' = ; 

it is obvious that one can treat the replacement as a unitary equivalence which does 
not alter spectral properties and preserves the system resonances (if there are any). 

7.2. Equivalence of resonance notions. In Scction[2]wc made it clear how 
important it is to establish connection between different objects labeled as reso- 
nances. Let us look now how this question looks like in the quantum graph setting. 
Let us begin with the resolvent resonances. One can write the resolvent of the 
graph Hamiltonian |PalOj . however, it is sufficient to inspect the spectral condi- 
tion encoded in it and its behavior in the complex plane. 

We employ an external complex scaling in which the external part are the 
semi- infinite leads where the functions are scaled as gje{x) = e^^'^gj{xe^) with an 
imaginary 9; as usual this rotates the essential spectrum of the transformed (non- 
selfadjoint) Hamiltonian into the lower complex halfplane and reveals the second- 



sheet poles. In particular the 'exterior' boundary values, to be inserted into (7.2 1, 
can be for gj{x) = Cj e'^'^^ written as gj{0) ~ eT^I'^gjs and 5j(0) — ike~^^'^gje 
with an appropriate gjg. On the other hand, the internal part of the graph is left 
unsealed. Choosing the solution on the j-th edge in the form fj{x) = aj sinfcx -I- 
bj cos kx we easily find its boundary values; for x = it is trivial, for x = Ij we 
use the standard transfer matrix. This allows us to express both {fj{0),fj{lj))'^ 
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and {fj{0),—fj{lj))'^ through the coefficients aj,bj, cf. [ELlOj and note the sign 



of fjilj) reflecting the fact that the derivatives entering (7.2 1 are outward ones. 



Inserting these boundary values into the coupling condition we arrive at the system 
/ ai \ / ai \ 



(7.3) iU-I)Ciik) 



ai 
bi 



519 



ik{U + I)C2{k) 



ai 
bi 

02 



5 AT 



where we have set Ci{k) = Ai&g {C^^\k), ' [k) 

diag (Cf ^(fc), ^(fc), . . . , d^\k),llMy.M) with 



,c[^Hk),lMxM) and C'2 = 



C['\k) = 





sin kL 



2 
1 

cos kL 



d^\k) = 



1 

cos kh 





sin kli 



and /, 



MxM being the M x M unit matrix. The solvability condition of the system 
(7.3) determines eigenvalues of the scaled non-selfadjoint operator, and mutatis 



mutandis, poles of the analytically continued resolvent of the original Hamiltonian. 

Looking at the same system from the scattering point of view we use the same 
solution as above on the internal edges while on the leads we take appropriate 
combinations of two planar waves, gj = Cje~'^^^ + dje^''^ . We look for the on-shell 
S-matrix S = S{k) which maps the vector of amplitudes of the incoming waves 
c = {c„} into the vector of amplitudes of the outgoing waves d — {dn}, and ask 
about its complex singularities, det — 0. This leads to the system 

f ai \ f ai \ 

bi 



(f/-/)Ci(fc) 



0-2 

bN 
ci + di 

\ CM + dn ) 



+ ik{U + I)C2{k) 



ai 
bi 
a2 



bN 
di - ci 

\ dM - CM J 



= 0; 



eliminating the variables aj, bj one can rewrite it as a system of M equations 
expressing the map S^^d = c. The condition under which the latter is not solvable, 
which is equivalent to our original question since S is unitary, reads 



(7.4) 



det [{U - I) Ci{k) + ik{U + I) C2{k)] =0, 



however, this is nothing else than the condition of solvability of the system (7.3 1. 
Thus we can make the following conclusion |EL10] . 

Theorem 7.1. The notions of resolvent and scattering resonances coincide for 
quantum graph Hamiltonians described above. 

Before proceeding further, let us mention one more way in which the resonance 
problem on a graph can be reformulated. To this purpose we rearrange the matrix 
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U permuting its rows and columns into the form U = (^^ [^^) , where Ui is the 
2N X 2N square matrix referring to the compact subgraph, U4 is the M x M square 
matrix related to the exterior part, and U2 and U3 are rectangular matrices of the 
size M X 2N and 2N x M, respectively, connecting the two. The system (7.3) can 
be then rewritten by elimination of the lead variables |EL10| as 

(7.5) {ij{k)~ I)F + i{U{k)+I)F' ^Q, 

where F := (/i, . . . , f2N)'^ : and similarly for F', are the internal boundary values, 
and the effective coupling matrix appearing in this condition is given by 

(7.6) U{k) = C/i - (1 - k)U2[{l - k)Ui - (fc + l)/]-i[/3 . 

In other words, we have been able to cast the problem into the form of spectral 
question for the compact core of the graph with the effective coupling replacing the 
leads by the non- unitary and energy-dependent matrix (7.6 1. 

7.3. Line with a stub. Next we will present several simple examples. In the 
first one F is a line to which a segment is attached at the point x — {). The Hilbert 
space is thus % := L^(M) © i^(0, 1) and we write its elements as columns -0 = 
To make the problem more interesting we suppose that the particle on the stub is 
exposed to a potential; the Hamiltonian acts 

{HijUx) = -f"{x) , {H^)2{x) - i-u" + Vu){x) , 

outside the junction, where V € L\^^{Q,l) having finite limits at both endpoints 
of the segment so — if the domain consists of functions vanishing in the vicinity 
of the junction — the corresponding deficiency indices are (3,3). The admissible 
Hamiltonians will be identified with self-adjoint extensions which, as before, can 
be conveniently characterized by appropriate boundary conditions. We will not 
explore all of them and restrict our attention to a subclass of those having the line 
component of the wave function continuous at the the junction, namely 

/(0+) = /(0-)=:/(0), «(0)-6/(0) + cu'(0), 

/'(0+) - /'(0-) = rf/(0) - W(0) , u{i) = ; 

at the free end of the stub we fix Dirichlet condition. The coefficient matrix K. = 
id -b) real; we restrict our attention to time-reversal invariant dynamics. The 
operator specified by the conditions (7.7) will be denoted as H/c- The parameter 
b controls the coupling; if & = the graph decomposes into the line with a point 
interaction at x = and the stub supporting Schrodinger operator he :— —353- + V 
with Robin condition u(0) = cu'(O) at the junction referring again to a; = 0. 

Let us begin with the scattering. To find the on-shell S-matrix we use the 
standard Ansatz, f{x) = e^^^ -\- r e"'*^^ and te^^^ on the line for x < and 2; > 0, 
respectively, while on the stub we take u{x) = j3ui{x) where ui is a solution to 
— m" + Vu = k^u corresponding to the boundary conditions ui(l) — 0, unique up to 
a multiplicative constant. Using the coupling conditions (|7.7|) we find 



2ikD{k) ' "^^'^ 2ikD{k) ' 

where 2ikD{k) :— 6^u^(0) + (d — 2ik){cu'^ — ui){Q)\ it is easy to check that these 
amplitudes satisfy |i(fc)p -I- |r(A:)p = 1. We note that H]c can have also an isolated 
eigenvalue; this happens if D{iK) — with k > 0. If 6 = such an eigenvalue exists 
provided d < and equals — jii^ ; it remains isolated for |&| small enough. 
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It is also not difficult to find the resolvent of Hjq. The tool is as usual Krcin's 
formula; we choose for comparison H)q^ corresponding to /Co = 0. In that case the 
operator decomposes, the kernel of line part being Gi(x,?/;z) = _L. g^^k-j/l -^vhere 
k := -y/z as usual; the stub part is — ""[^^^"m,'^)^'' : where ui has been introduced 
above, uq is similarly a solution corresponding to the condition uo(0) ~ cuq(O), and 
W{uq,ui) is the Wronskian of the two functions. The sought kernel then equals 

(Hfc ~ z)-\x,y) = {Hfc„ - z)-\x,y) + ^ X,kik)F,{x)Fk{y) , 

3 = 1,2 

where the vectors Fj can be chosen as Fi{x) := and ^2(2^) '■— („,(a;)); 

note that the stub part vanishes at a; = 0. The coefficients are obtained from the 
requirement that the resolvent must m ap any vector of H into the domain of Hjc] 
a straightforward computation (ES94j gives 



together with A12 (A:) = A2i(A:) = hD{k)-^. We see, in particular, that the coefficient 
denominator zeros in the complex plane coincide with those of the on-shell S- matrix 
as we expect based on Theorem 7.1 proved above. 

In the decoupled case, 6 = 0, the expression for D{k) factorizes giving rise 
to eigenvalues of the operator he introduced above which are embedded in the 
continuous spectrum of the line Hamiltonian; the coupling turns them generally 
into resonances. In the case case of weak coupling, i.e. for small |6| one c an solv e 
the condition D{k) = perturbatively arriving at the following conclusion ES94j . 

Proposition 7.2. Let fc„ refer to the n-th eigenvalue of he and denote by Xn 
the corresponding normalized eigenf unction; then for all sufficiently small \b\ there 
is a unique resolvent pole in the vicinity of fc„ given by 

This gives, in particular, the inverse value of resonance lifetime in the weak- 
coupling case, Imz„(6) = i6^Xn(0)^(2fcn + id)^^ + 0(6"*). The simple form of the 
condition D{k) — allows us, however, to go beyond the weak coupling and to 
trace numerically the pole traject ories as the coupling constant b runs over the 
reals. Examples are worked out in ES94j but we will not describe them here and 
limit ourselves with mentioning an important particular situation. 

It concerns the case when the motion in the stub is free and the decoupled 
operator is specified by the Dirichlet boundary condition, i.e. V = Q and c = d = 0. 
The condition D{k) ~ can be then solved analytically. Indeed, the embedded 
eigenvalues are fc^ with kn ^ and the equation reduces to tanH — — ^ solved 

by 



^ + ^lny, ... \b\<V2 



(7.8) kn{b) 

Hence the poles move with the increasing |6| vertically down in the fc-plane and for 
|6| > y/2 they ascend, again vertically, returning to eigenvalues of Neumann version 
of he as |5| — >■ 00. An important conclusion from this example is that poles may 
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disappear to infinite distance from tlie reai axis and a quantum graph may have no 
resonances at all, as it happens here for |&| = \/2- 

7.4. Regeneration in decay: a lasso graph. Let us next describe another 
simple example, now with a lasso-shaped F consisting of a circular loop of radius R 
to which a halfline lead is attached. This time we shall suppose that the particle is 
charged and the graph is placed into a homogeneous magnetic field of intensity B 
perpendicular to the loop plan^ The vector potential can be then chosen tangent 
to the loop with the modulus A = ^ BR = * , where $ is the flux through the 
loop and L is its perimeter. With the convention we use, e = c = 2m = h = 1, the 
natural flux unit is ^ = 27r, so we can also write A = ^ where (p is the flux value in 
these units. The Hilbert space of the lasso-graph model is H := L'^{0, L) ® 
the wave functions are written as columns, ijj = (^) . 

To construct the Hamiltonian we begin with the operator describing the free 
motion on the loop and the lead under the assumption that the graph vertex is 
'fully disconnected', in other words Hoo = ^fioop(-B) ffl ^^haiflinc, where 

/ . d \2 d^ 

Hloop{B) ^ \ * + "^j ' ^halfiino = ~ 

with Dirichlet condition, u{0) = u{L) = /(O) = at the junction. The spectrum 
of Hioop is discrete of multiplicity two; the eigenfunctions Xn{x) = sin ( jg) 

with n — 1, 2, . . . correspond to the eigenvalues (^)^ which are embedded into the 
continuous spectrum of iJhaifiinc covering the interval [0, oo); note that the effect 
of the magnetic field on the disconnected loop amounts to a unitary equivalence, 
Hyoop{B) = U^AHioopiO)UA where := e*^-u(a;). 

Restricting the domain of to functions vanishing in the vicinity of the 
junction we get a symmetric operator with deficiency indices (3,3). We are going 



to consider a subclass of its self-adjoint extensions analogous to (7.7 1 characterized 
by three real parameters; the Hamiltonian will act as 

on functions from W^^^{T) continuous on the loop, u(0) = u{L), which satisfy 

(7.9) f{0)^LuuiO)+^if'{0), u'{0)-u'iL)^au{0)-Ljf'{0), 

for some Q;,/i,a; G M the latter being the coupling constant. This includes a partic- 
ular case of S-coupling corresponding to ^ = and w = 1 in which case the wave 
functions are fully continuous, 

(7.10) u{0)^u{L)^f{0), u'(0)-u'(L)+/'(0) =«/(0); 

in the fully decoupled case we have a = cx) as the notation suggests. For simplicity 
we will write -ffa^o.i = Ha- Note that in general the vector potential enters the 
coupling conditions |KS03j but here the outward tangent components of A at the 
junction have opposite signs so their contributions cancel mutually. 

Let us start again with scattering, i.e. the reflection of the particle traveling 
along the halfline from the magnetic-loop end. To find the generalized eigenvectors. 



■^The assumptions of homogeneity and field direction are here for simplicity only, in fact the 
only thing which matters in the model is the magnetic flux through the loop. 
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Ha,n,u}{B)ip ~ fc^Vj we use u(x) — f5e ''^^ sin(/cx + 7) and f{x) — e -|_ 7- e*'=^ as 



the Ansatz; using the coupUng conditions (7.9 1 we get after a simple algebra 



^^^^ ^ _ (1 + ik^i) [g - ^ (cos $ ~ cos kL) ] + lu^k 
(1 — ik^) [a — ^^j^ (cos $ — cos kL) ] — iw^fc 

for the reflection amplitude. The Hamiltonian Ha,fi,ui{B) can have also isolated 
eigenvalues but we shall skip this effect referring to |Ex97| . On the other hand, 
it is important to mention that there may exist positive eigenvalues embedded in 



the continuous spectrum even if a; 7^ 0. In view of (7.9) it is possible if u(0) = 
u'{0)— u'{L) = 0, hence such bound states exist only at integer/halfinteger values 
of the magnetic flux (in the natural units) and the corresponding eigenfunctions 
are the Xn's mentioned above with even n for integer and odd n for (p halflnteger. 

Next we find the resolvent of H^.^i^uiiB) using again Krein's formula to compare 
it to that of with the kernel diag (e-^^f--^) ""^"</2t^^>-^\ sinfc.< c"'-> 
The sought resolvent kernel can be then written as 

2 

{x,y;k) ^ Goo{x,y;k) + ^ \.ji{k)Fj{x)Fi{y) , 
where the deficiency subspaces involved are chosen in the form 



„ , , I w(x)\ „ , X / 



ikx 



with w{x) \— e*^^ ^^"fin/eL^^^^ coefficients \ji{k) given by |Ex97| 

1 - 7 ,, To/- cos$ — cosfcL ^1 , ,2 



•^11 ^ „ N , A22 — 



D{k) ' D{k) 
together with A12 = A21 = ^jjp^i where 

cos $ — cos kL 



D{k) = D{a, /iw; k) := (1 — i/ifc) 



2k- 



sin kL 



iuj k . 



In the case of a 5-coupling, in particular, the coefficients acquire a simple form, 
Xji{k) = —D{k)~^, j,l = 1,2 . As expected, the denominator D{k) determining the 
singularities is the same as for the on-shell S-matrix. A simple form of the condition 
D{k) = allows us to follow the pole trajectories, both with respect to the coupling 
parameters and the flux <&. At the same time, knowing the resolvent of Ha^^^u{B) 
we can express the decay law for states supported at the initial moment t = on 
the loop only; we will not go into details and refer the reader to |Ex97| where the 
appropriate formulae and graphs are worked out. 

Let us just mention one amusing feature of this model which can be regarded 
as an analogue of the effect known in particle physics as regeneration in decay of 
neutral kaons and illustrates that intuition may misguide you when dealing with 
quantum systems. Consider the lasso graph with the initial wave function u on 
the loop such that x 1— )• e*"^^u(a;) has no definite symmetry with respect to the 
connection point a; = 0. If the flux value (/> is integer, the A-even component 
represents a superposition of embedded-eigenvalue bound states mentioned above, 
thus it survives, while the A-odd one dies out. Suppose that after a sufficiently long 
time we decouple the lead and attach it at a different point (or we may have a loop 
with two leads which may be switched on and off independently). For the decay 
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of the surviving state the symmetry with respect to the new junction is important; 
from this point of view it is again a superposition of an A-even and an A-odd part, 
possibly even with same weights if the distance between the two junctions is jL. 

7.5. Resonances from rationality violation. The above simple examples 
illustrated that resonances are a frequent phenomenon in quantum graph models. 
To underline this point we shall describe in this section another mechanism giving 
rise to resonances, this time without need to change the coupling parameters. The 
observation behind this claim is that quantum-graph Hamiltonians may have em- 
bedded eigenvalues even if no edges are disconnected which is related to the fact 
that the unique continuation principle is generally not valid here and one can have 
compactly supported eigenfunctions. Indeed, eigenfunctions of a graph Laplacian 
are trigonometric functions, hence it may happen that the graph has a loop and the 
vertices on it have rationally related distances such that the eigenfunction has zeros 
there and the rest of the graph 'does not know' about it. Let us present briefly two 
such examples referring to |EL10j for more details. 

7.5.1. A loop with two leads. In this case F consists of two internal edges of 
lengths ^1, I2 and one halfline attached at each of their endpoints, corresponding to 
the Hilbert space is L'^{R+) (B L^(^'^) O L'^{[0,li]) (B L'^{[0,l2]); states of the system 
are correspondingly described by columns ip = {gi, g2, fi, f2)'^ ■ The Hamiltonian 
H is supposed to act as negative Laplacian, ^ 1-^ — V'"; separately on each edge. 



We consider coupling analogous to (7.7) assuming that the functions of DomiJ are 
continuous on the loop, /i(0) = /2(0) and fi{li) = /2(^2), and satisfy 

/l(0) = a^\fm + + 713; (0) , fl{h) = -0^2\flih) + /2(^2)) + 72.g2(0) , 

31 (0) = 7l(/{(0) + m) + «r'3i(0) , 52(0) = -72(/{(^l) + /2(^2)) + «2"'<?2(0) , 

for some a^, Uj, jj G M. Since we want to examine behavior of the model with 
respect to the lengths of internal edges, let us parametrize them as h = /(I — A), I2 = 
1{1 + A) with A e [0,1]; changing A thus effectively means moving one of the 
connections points around the loop from the antipolar position for A = to merging 
of the two vertices for A = 1. Due to the presence of the semi-infinite leads the 
essential (continuous) spectrum of H is [0, 00). If we consider the loop itself, it 
has a discrete spectrum consisting of eigenvalues fc^ where kn — ^ with n G Z. 
The corresponding eigenfunctions have nodes spaced by - for n 7^ 0, hence H has 
embedded eigenvalues if the leads are attached to the loop at some of them. 

If the rationality of the junction distances is violated these eigenvalues turn 
into resonances. The condition determining the singularities can be found in the 
same way as in the previous examples; it reads 



4fc' 

sin fc;(l - A) sin kl{l + A) - sin"' kl + k 



Pi{k)h{k) lpi{k) P2{k) 



1 1 



sin2fcZ = 0, 



where l3-\k) a'^ + ^Z^- One can solve it perturbatively but also to find 

numerically its solution describing pole trajectories as A runs through [0,1]. The 
analysis presented in |EL10j shows that various situations may occur, for instance, 
a pole returning to the real axis after one or more loops in the complex plane — 
an example is shown in Fig. [3] — or a trajectory ending up in the lower halfplane 
at the endpoint of the parameter interval. 
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Figure 3. Trajectory of the resonance pole in the momentum 
plane starting from fcg = 27r corresponding to A = for / = 1 
and the coefficients values a^^ — 1, a^^ — —2, ji = 1, Oj^^ = 0, 
(5^^ = 1, 72 = 1, n = 2. 



7.5.2. A cross-shaped graph. We add one more simple example to illustrate that 
the same effect may occur even if the graph has no loops. Consider a cross-shaped 
r consisting of two leads and two internal edges attached to the leads at one point; 
the lengths of the internal edges will be li = 1{1 — X) and I2 = 1{1 + X). The 
Hamiltonian acts again as — d^/dcc^ on the corresponding Hilbert space L^(]R+) © 
I/2(R+) © L^([0,li]) © L^([0, Z2]) the elements of which are described by columns 

= (5ii52, fi, f2)^- For simplicity we restrict our attention to the 5 coupling at the 
vertex and Dirichlet conditions at the loose ends, i.e. /i(0) — /2(0) = 5i(0) = (?2(0) 
and fi{li) — f2{h) — together with the requirement 

/{(0) + /2(0)+.9i(0) + <?2(0)=«/i(0) 

for a G M. In the same way as above we can derive resonance condition in the form 
k sin 2kl + {a — 2ik) sin kl{l — A) sin fcZ(l + A) = 0, or equivalently 

2k sin 2kl + (a - 2ik) (cos 2klX - cos 2kl) = . 

Asking when the solution is real we note that this happens if the real and imaginary 
parts of the left-hand side vanish. We find easily that it is the case if A = 1 — 2m /n , 
No 9 m < n/2, while if this rationality relation is violated the poles move into 
the lower halfplane. The condition can be again solved numerically giving pole 
trajectories for various parameter values; in addition to the possibilities mentioned 
above we can have trajectories returning to different embedded eigenvalues — an 
example shown in Fig. |4] calls to mind the effect of quantum anholonomy |Ch98| 
— as well as avoided trajectory crossings, etc., see [EL 101 for more details. 

7.5.3. Local multiplicity preservation. Let us turn from examples to the gen- 
eral case and consider an eigenvalue kg with multiplicity d of a quantum-graph 
Hamiltonian H which is embedded in the continuous spectrum due to rationality 
relations between the edges of F. We consider graphs F^ with modified edge lengths 
I'j = loirij + Ej) assuming that € N for j G {1, . . . , n} while Uj may not be an 
integer for j G {n -I- 1, . . . , A^} where N := The analysis of the perturbation is 
a bit involved, see [ELlOj for details, leading to the following conclusion. 
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Figure 4. Resonance pole trajectory for a = 1 and n — 2. 



Theorem 7.3. Let T be a quantum graph with N finite edges of the lengths 
Ij, M infinite edges, and the coupling described by the matrix U — (^^ ^|), where 
Ui corresponds to the coupling between the infinite edges. Let fco > correspond 
to a pole of the resolvent (H — /cq)~^ of multiplicity d. Let he a geometrically 
perturbed quantum graph with edge lengths loijij +£j) described above and the same 
coupling as V . Then there exists an Eq > such that for all e £ ^^^(O) the sum of 
multiplicities of the resolvent poles in the vicinity of ko is d. 

8. High-energy behavior of quantum-graph resonances 

Now we will look at quantum-graph resonances from a different point of view 
and ask about asymptotics of their numbers at high energies. Following the papers 
[DPlll IDELIO) we are going to show, in particular, that it may often happen 
that this asymptotics does not follow the usual Weyl's law. Following the standard 
convention we will count in this section embedded eigenvalues among resonances 
speaking about the poles in the open lower halfplane as of 'true' resonances. 

8.1. Wey l asymptotics criterion. It is useful for our purpose to rewrite the 
condition (7.5) in terms of the exponentials e**^'^ and e"*'^'^ using for brevity the 



symbols ef :— o^^^h and e^ :— Ifj^^^e^ — e^^''^ , where V = YTj=i h is the size of 
the finite part of F. The condition then becomes 

F{k) := det |i[(C/-/) + k{U+I)]E^{k) + \[{U-I) - k{U+L)]E2{k) 

+k{U+I)E3 + {U-I)Ei + [([/-/) - k{U+I)] diag (0, . . . , 0, /mxm) | = , 

where E.,{k) diag (^E^^^\ e'i^\ . . . , £;|^\ 0, . . . , , i = 1, 2, 3, 4, are matrices con- 
sisting of N nontrivial 2x2 blocks 

and a trivial M x M part. To analyze the asymptotics we employ the usual counting 
function N(R,F) defined for an entire function F{-) by 

N{R, F) = tt{fc : F{k) = and |fc| < R} , 
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where the algebraic multiplicities of the zeros are taken into account. With the 
above spectral condition in mind we have to apply it to trigonometric polynomials 
with rational-function coefficients. We need the following result [DELIO] which is 
a simple consequence of a classical theorem by Langer |La31) . 

Theorem 8.1. LetF{k) =J2r=o'^r{k)e'''''- , wh ere ar{k) are rational functions 
oj the complex variable k with complex coefficients, and ct^ G I^, < < • • • < (T„. 
Suppose that Yimk^ao o,o{k) ^ and limfc_>.oo an(fc) 7^ 0. There exists a compact set 
C C, real numbers and positive K^, r — 1, . . . ,n, such that the zeros of F{k) 
outside Q lie in one of n logarithmic strips, each one bounded between the curves 
— Imfc + rUr log = ±Kr- The counting function behaves in the limit i? — > cx) as 

N{R,F) = + 
vr 

To apply this result it is useful to pass to effective energy-dependent coupling 



(7.6) which makes it possible to cast the spectral condition into a simpler form, 



(8.1) F{k) ^ dct S^^[{U{k) - I) + k{U{k) + T)]Ei{k) 

+ l[{U{k) -I)- k{U(k) + I)]E2{k) -f k{i}{k) + 7)^3 + {U{k) - /)^4| - , 

where Ej are the nontrivial 2 x 27V parts of the matrices Ej , the first two of them 
being energy-dependent, and / denotes the 27V x 27V unit matrix. Then we have 
the following criterion I DELlOj for the asymptotics to be of Weyl type. 

Theorem 8.2. Assume a quantum graph (T,Hu) corresponding to T with 
finitely many edges and the coupling at vertices Xj given by unitary matrices Uj. 
The asymptotics of the resonance counting function as i? — > oo is of the form 

2W 

N{R,F) = R + 0{1), 

TT 

where the effective size of the graph W satisfies < W < V := X^jLi h- More- 
over, W < V holds if and only if there exists a vertex where the corresponding 
energy- dependent coupling matrix Uj{k) has an eigenvalue or for all k. 

To prove the theorem one has to realize that cr„ — V and ctq — —V, hence 
the asymptotics is not of Weyl type iff either the senior or the junior coefficient 
in expression of F{k), i.e. those of e^, vanish. By a straightforward computation 

[DELlOj we find that they equal (f)^ det [iU{k) -I)± k{U{k) + /)], respectively, 
and therefore they vanish under the condition stated in the theorem. 

Before proceeding further, let us mention that the asymptotic number of res- 
onances is not the only thing of interest. One can investigate other asymptotic 
properties such as the distribution of resonance pole spacings; quantum graphs are 
known to be a suitable laboratory to study quantum chaotic effects |KoS03| . 

8.2. Permutation-symmetric coupling. Let us first look what the above 
criterion means in a particular class of vertex couplings which are invariant with 
respect to permutations of the edges connected at the vertex. It is easy to see that 
such couplings are described by matrices of the form Uj = ajJ-\- bjl, where a-,-, 
bj are complex numbers satisfying \bj\ = 1 and \bj -\- ajdegXj\ = 1; the symbol J 
denotes the square matrix all of whose entries equal to one and I stands for the 
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unit matrix. Important examples arc the 5-coupling analogous to (7.101, with the 
functions continuous at the vertex and the sum of outward derivatives proportional 
to their common value, corresponding to Uj = 'j~^j^J~ I, where dj is the number 
of edges emanating from the vertex Xj and aj e M is the coupling strength, and the 
SI. -coupling corresponding to Uj — — ^ J + / with Pj e M for which the roles of 
functions and derivatives are interchanged. The particular cases aj = and /Sj = 
are usually referred to as the Kirchhoff and anti-Kirchhoff condition, respectively. 

Consider a vertex which connects p internal and q external edges. For matrices 
of the form Uj = ajJ + bjl it is an easy exercise to invert them and to find the 
effective energy-dependent coupling; this allows us to make the following claim. 

Theorem 8.3. Let (T,Hij) be a quantum graph with permutation- symmetric 
coupling conditions at the vertices, Uj = ajj + bji. Then it has a non-Weyl asymp- 
totics if and only if at least one of its vertices is balanced in the sense that p ^ q, 
and the coupling at this vertex satisfies one the following conditions: 

(a) fm = /„, ym,n < 2p, /m = 0. i.e. U ^ ~ h. 



(b) fL ^fn, ym,n< 2p, E™=i fj =0, ».e. f/ = J2px2p + h 



'2px2p ; 
'Jj'J2py.2p -|- I2px2p ■ 

In other words, if the graph has a balanced vertex there are exactly two situations 
when the asymptotics is non-Weyl, either if the coupling is Kirchhoff — which is 
the case where the effect was first noted in |DP11) — or if it is anti-Kirchhoff. 

8.2.1. An example: a loop with two leads. To illustrate the above claim let us 
return to the graph of Example |7.5.1[ It is balanced if the two leads are attached at 
the same point. Changing slightly the notation we suppose that the loop length is I 
and consider negative Laplacian on the Hilbcrt space is L^(0, /) ® ® L^(M+) 

with its elements written as (w, /i, f2)^ defined on functions from Wi^^(T) satisfying 
the requirements u(0) = /i(0) and u{l) = /2(0) together with 

au{Q) = u'{0) + f[ (0) l3{-u'{l) + /^(O)) , au{l) = /3{u'{0) + /( (0)) - u\l) + /^(O) 

with real parameters a and /3; the choice (3 — 1 corresponds to the 'overall' 6- 
coupling of strength a, while (5 — decouples two 'inner-outer' pairs of mutually 
meeting edges turning F into a line with two ^-interactions at the distance I. In 
terms of the quatities = e**'^' the pole condition can be written ^DELlOj as 

ia'^e+ -\- 4:ka(3 - i[a{a - Aik) -f 4:P{(3'^ - 1)] e" _ 
4(/32 - l) + a{a-4i) ^ 

The coefficient of 6+ vanishes iff a — 0, the one in the second term for /3 = or if 
\(3\ 7^ 1 and a = 0, while the coefficient does not vanish for any combination of 
a and /3. The graph has thus a non-Weyl asymptotics iff a — 0. If, in addition, 
1/3 1 ^ 1, then all resonances are confined to a circle, i.e. the graph has zero 'effective 
size'. The only exceptions are the Kirchhoff condition, /3 = 1 and a = 0, and its 
anti-Kirhhoff counterpart, (i — —1 and a = 0, for which one half of the resonances 
is asymptotically preserved, in other words, the effective size of the graph is ^l. 

We can demonstrate how a 'half of the resonances disappears using the example 
of the (5-coupling, (3—1. The resonance equation in this case becomes 

— asinfcZ -I- 2fc(l -I- isinfc/ — cosfc^) ^ 

a — 4i 

A simple calculation shows that the graph Hamiltonian has a sequence of embedded 
eigenvalues with k — n e Z, and a family of resonances given by solutions 
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Figure 5. Graph with a balanced vertex 



to e = — 1 + The former do not depend on a, while the latter behave like 

Imfc = -y In- +0(1) , Rcfc = n7r + e'(a), 
I a 

as a — > 0, hence all the 'true' resonances escape to the imaginary infinity in the 



limit, in analogy with the similar pole behavior described by relation (7.8) 



8.3. The mechanism behind a non-Weyl asymptotics. One naturally 
asks about reasons why graphs with balanced vertices and Kirchhoff/anti-Kirchhoff 
coupling have smaller than expected effective size. A simple observation is that if 
such a vertex has degree one, then KirchhofF coupling between an external and 
internal edge is in fact no coupling at all, hence the internal edge can be regarded 
as a part of the lead and the effective size is diminished by its length. We are going 
to show that this remains true in a sense also when the degree is larger than one. 

8.3.1. Kirchhojf 'size reduction'. A graph with a balanced vertex can be always 
thought as having the shape sketched in Fig. [5] with a vertex Xi which connects 
p internal edges of the same length Iq and p external edges; if the internal ones 
have different lengths we simply add a needed number of 'dummy' vertices. We 
will suppose that the coupling at Xi is invariant with respect to edge permutations 
being described by a unitary matrix U^^^ = aJ2px2p + bl2px2p', the coupling of the 
other internal edge ends to the rest of the graph, denoted here as Fq, is described 
hy SL q X q matrix U^'^\ where q > p (which may express also the topology of Fq). 

To find how the effective size of such a quantum graph may look like we employ 



the following property which can be derived easily from coupling condition (7.2) 



Proposition 8.4. LetT be the graph described above with the coupling given by 
arbitrary U^^^ andU^^\ Let further V be an arbitrary unitary pxp matrix, V^^"^ := 
diag (y, y) and F'-^' :— diag (/(g_p)x(g-p)i 1^) be 2p x 2p and qx q block diagonal 
matrices, respectively. Then H on T is unitarily equivalent to the Hamiltonian 
Hy on the graph with the same topology and the coupling given by the matrices 
and [F(2)j-ij/(2)y(2)^ respectively. 

Application to the couplings described by U^^^ — aJ2px2p + bl2px2p at Xi is 
straightforward. One has to choose the columns of V as an orthonormal set of 
eigenvectors of the corresponding pxp block aJpxp + bipxp of U^-^\ the first one of 
them being ^ (1, 1, . . . , 1)"^. The transformed matrix decouples 
then into blocks connecting only the pairs {vj,gj). The first one of these, corre- 
sponding to a symmetrization of all the Uj's and /^'s, leads to the 2x2 matrix 
C^2x2 = 0,pJ2x2+bl2x2, whilc the other lead to separation of the corresponding inter- 
nal and external edges described by Robin conditions {b~l)vj{0) + i{b+l)vj{0) = 
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Figure 6. A polygonal balanced graph. 



and (6 — l)5j(0) + i{b + 1)5^(0) = for j = 2, . . . ,p. We note that it resembles the 
reduction procedure of a tree graph due to Solomyak jSS02| . 

It is easy to see that the 'overall' Kirchhoff/anti-KirchhofF condition at Xi 
is transformed to the 'line' KirchhofF/anti-Kirchhoff condition in the subspace of 
permutation-symmetric functions, leading to reduction of the graph effective size as 
mentioned above. In all the other cases the point interaction corresponding to the 
matrix apJ2x2 + bl2x2 is nontrivial, and consequently, the graph size is preserved. 

8.3.2. Global character of non-Weyl asymptotics. The above reasoning might 
lead one to the conclusion that the effect discussed here is of a local character. We 
want to show now that while this is true concerning the occurrence of non-Weyl 
asymptotics, the effective size of a non-Weyl quantum graph is a global property 
because it may depend on the graph F as a whole. 

We will use an example to justify this claim. We shall consider the graph F„ 
with an integer n > 3 which contains a regular n-gon, each edge of which has length 
I. To each of its vertices two semi- infinite leads are attached, cf. Fig. [6] Hence all 
the vertices of F„ are balanced, and if the coupling in them is of Kirchhoff type the 
effective size Wn of the graph is by Theorem |8.3| strictly less than the actual size 
Vn = n£. It appears, however, that the actual value of the effective size depends in 
this case on the number n of polygon vertices. 

Since all the internal edges have the same length, the system has a rotational 
symmetry. One can thus perform a 'discrete Floquet' analysis and investigate cells 
consisting of two internal and two external edges; the wave functions at the ends 
of the former have to differ by a multiplicative factor uj such that uj" — 1. After a 
simple computation |DEL10| we conclude that there is a resonance at iff 



(8.2) 



I) 



-ikl 



0. 



The 'Floquet component' of H has thus effective size — if + 1 7^ 
while for + 1 = we have no resonances, W^^ — 0. Summing finally over all the 
Lj with = 1 we arrive at the following conclusion. 



Theorem 8.5. The effective size of the graph r„ with Kirhhoff coupling is 

^nl if n 7^ mod 4 

^{n-2)l if n = mod 4 
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Note that if one puts u = e^^ in (8.2) the resonance poles behave according to 



i(^ln(c 



where n S Z is arbitrary, hence they escape to imaginary infinity as 9 — > i^Tr. 
Of course, the Floquet variable is discrete, 9 = j = 0, 1, . . . ,n— 1, neverthe- 
less, the limit still illustrates the mechanism of the resonance 'disappearance'; it is 
illustrative to look at the behavior of the solutions for large values of n. 

8.4. Non-Weyl graphs with non-balanced vertices. Now we are going 
to show that there are many more graphs with non-Weyl asymptotics once we 
abandon the assumption of permutation symmetry of the vertex couplings. For the 
sake of brevity, we limit ourselves again to a simple example. In order to formulate 
it, however, we state first a general property of the type of Proposition |8.4| above. 
Specifically, we will ask what happens if the coupling matrix [/ of a guantum graph 
is replaced by W~^UW, where is a block diagonal matrix of the form 

w = I ^"^^pxp ^ 

' Wi 



and W4 is a unitary q x q matrix. The following claim is obtained easily from ( 7.2 ) 



Proposition 8.6. The family of resonances of a quantum- graph Hamiltonian 
Hjj does not change if the original coupling matrix U is replaced by W^^UW. 

Let us turn now to the example which concerns the graph investigated in 
Sec. |7.3[ a line with a stub of length I, this time without a potential. Changing 
slightly the notation we use the symbols fj for wave function on the two halflines 
and u for the stub. The function from the domain of any Hamiltonian Hu are 
locally W^'^ and satisfy the conditions u{l) + cu'{l) = with c e M U {00} and 

({/ - /) (^.(0), /i(0), /2(0))T + + /) {u\Q)J'MJm? = . 



We split r into two parts in a way different from Sec. |7.3| choosing the coupling 
described by C/q := diag(('|' p), e*''^) which gives two halflines with the conditions 
u{l) + cu'{l) — and /2(0) + cot^/2(0) = 0, respectively, at their endpoints; 
the flrst part consists of the halfline number one and the stub joined by Kirchhoff 
coupling. It is obvious that such a graph has at most two resonances, and thus a 
non-Weyl asymptotics. We now replace C/q by Uw = W~^UqW with 





\/l — e'''^^ „^(5j(v2+¥'3-¥'i) 

for some r G [0, 1] and obtain for every fixed value of ip and c a three-parameter 
family of coupling conditions described by the unitary matrix 

/ re^'^i Vl ~ r'^e'f^ 

Uw = re-"^' (1 - r2)e*'^ -r^l r'^e-''^-''f'+'f^~'P^^ 

each of which has the same resonances as Uq by Proposition |8.6[ The associated 
quantum graphs are thus of non-Weyl type despite the fact that their edges are 
connected in a single vertex of F which is not balanced. 
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Note that among these couphngs one can find, in particular, the one mentioned 



above in connection with relation (7.8); choosing ip = tt and c — 0, and furthermore, 
(pi — '~P2 = and r = 2^^/^, we get the conditions 

/i(0) = /2(0), u{0)^V2h{0): fm-fm^~V2u'{0), 



or (7.7 1 with b = \/2 and c = d = 0. Similarly, such conditions with b — —\f2 and 
c — d — correspond lo Lp\ — — and r = 2~^/^; both these quantum graphs 
have no resonances at all. This fact is easily understandable, for instance, if we 
regard the line with the stub as a tree with the root at the end of the stub and 
apply the Solomyak reduction procedure [SS02 mentioned above. 

8.5. Magnetic field infiuence. Let us finally look how can the high-energy 
asymptotics be influenced by a magnetic field. We have encountered magnetic 
quantum graphs already in the example of Sec. |7.4[ now we look at them in more 
generality. We consider a graph F with a set of vertices and set of edges 
containing N finite edges and M infinite leads. We assume that it is equipped with 
the operator H acting as on the infinite leads and as — (g^ ^iAj(x)Y on the 

internal edges, where Aj is the tangent component of the vector potential; without 
loss of generality we may neglect it on external leads because one can always remove 
it there by a gauge transformation. The Hamiltonian domain consists of functions 
from Wj^f (r) which satisfy (C/, -7)^'^ + /)(^'^- = at the vertex Xj. 

As before it is useful to pass to a graph Fq with a single vertex of degree {^N + M) 
in which the coupling is described by the condition 

([/ - /)* + i{U + /)(*' + iA^) = , 

where the matrix U consists of the blocks Uj corresponding to the vertices of F 
and the matrix A is composed of tangent components of the vector potential at the 
vertices, A = diag (Ai(0), -Ai{h),..., ^a,(0), 0, . . . , 0). 

Using the local gauge transformation Tpj{x) 'ijjj{x)e~^^^''^'' with Xji^)' = 
Aj{x) one can get rid of the explicit dependence of coupling conditions on the 

,2 

magnetic field and arrive thus at the Hamiltonian acting as — with the coupling 
conditions given by a transformed unitary matrix, 

(8.3) {UA-i)'^ + i{UA + iW = ih Ua--=:fuf-\ 

with J- = diag (1, exp {i'^i), . . . , 1, exp (j<i>jv), 1, . . . , 1) containing magnetic fluxes 
$j — Aj{x)dx. Furthermore, one can reduce the analysis to investigation of 
the compact core of F with an effective energy-dep end ent coupling described by the 



matrix ?7^(/c) obtained from Ua in analogy with (7.5). 



To answer the question mentioned above we employ another property of the 



type of Proposition 8.4 This time we consider replacement of U by V~^UV where 
V — Y^) is unitary block-diagonal matrix consisting of a 2N x 2N block Vi 
and an M X M block V2; resonances are again invariant under this transformation. 
With respect to the relation between U and Ua we get the following result [ELllj . 

Theorem 8.7. A quantum graph with a magnetic field described by a vector 
potential A is of non-Weyl type if and only if the same is true for A = 0. 

In other words, magnetic field alone cannot switch a graph with non-Weyl 
asymptotics into Weyl type and vice versa. On the other hand, the magnetic field 
can change the effective size of a non-Weyl graph. To illustrate this claim, let 
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US return to the example discussed in Sec. |7.5.1[ a loop with two external leads 
KirchhofF-coupled to a single point, now we add a magnetic field. It is straightfor- 
ward to check |EL11] that the condition determining the resonance pole becomes 

-2 cos $ + 6"''=' =0, 

where $ is the magnetic flux through the loop. The graph is non-Weyl as the term 
with e**^' is missing on the left-hand side; if $ = ±7r/2 (modTr), that is, for odd 
multiples of a quarter of the flux quantum 27r, the ^-independent term disappears 
and the effective size of the graph becomes zero. 

The conclusions of the example can be generalized |ELllj to any graph with 
a single internal edge: if the elements of the effective 2x2 coupling matrix satisfy 
|'Wi2(A;)| = |M2i(fc)| for any A; > there is a magnetic field such that the graph under 
its influence has at most finite number of resonances. 

9. Leaky graphs: a caricature of quantum v^rire and dots 

A different class of quantum graph models employs Schrodinger operators which 
can be formally written as —A — a6{x — T) where T C M.'^ is a graph; one usually 
speaks about 'leaky' graphs. Their advantage is that they can take into account 
tunneling between different parts of the graph as well as its geometry beyond just 
the edge lengths. A survey of results concerning these models can be found in 
[Ex08] . In particular, even a simple F like an infinite non-straight curve can give 
rise to resonances jEN03| . however, one needs a numerical analysis to reveal them. 

9.1. The model. Instead we will describe here a simple model of this type 
which can be regarded as a caricature description of a system consisting of a quan- 
tum wire and one or several quantum dots. The state Hilbert space of the model 
is and the Hamiltonian can be formally written as 

n 

-A - aS{x + - y^"'') , 

where a > 0, S := {(xi,0); xi G R"^}, and H := {y^^^jf^i C \ E. The formal 
coupling constants of the two-dimensional S potentials are marked by tildes to 
stress they are not identical with the proper coupling parameters /3i which we shall 
introduce below. Following the standard prescription |AGHH| one can define the 
operator rigorously |EK04j by introducing appropriated boundary conditions on 
sun. Consider functions V' & W^^^ {R'^\{T,UU))nL'^ continuous on E. For a smaU 
enough p > the restriction ijj fc^ ; to the circle Cp^i :— {q £ : \q ~ ?;'^*'| = p} is 
weU defined; we say that ip belongs to D{Ha,p) iff {dl^ + d^J^p on \ (E U H) 
belongs to in the sense of distributions and the limits 

:= - lirn ^^-0 tcp ^ , := lim[-0 \c,i +Si(?/')lnp] , z = 1, . . . ,n, 

p->-olnp p-i-o 

5s(V')(a;i) := dx2i'{xi,0+) - dj;^'ip{xi,0-) , ns{ip){xi) := i'{xi,0) 
exist, they are finite, and satisfy the relations 



(9.1) 2^As,(i^) = n^iip) , Ss(0)(xi) = 



-aV,^{'tp){xi) , 
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where /3i G M are the true couphng parameters; we put j3 = . . . ,/3„) in the 
foUowmg. On this domain we define the operator Ha.ji '■ D{Ha.i3) ^ i^(M^) by 

Ha^,3^{x) = -Aip{x) for 2;eM2\(I]un). 

It is now a standard thing to check that Ha,i3 is essentiahy self-adjoint |EK04) : we 
shah regard in the following its closure denoted as Ha^p as the rigorous counterpart 
to the above mentioned formal model Hamiltonian. 

To find the resolvent of i?a,/3 we start from R{z) = (—A — z)^^ which is for any 
z G C \ [0, oo) an integral operator with the kernel Gz(x, x') — j^Ko{y/~z\x — x'\), 
where Kq is the Macdonald function and z i— ^ ^/z has conventionally a cut along 
the positive halfline; we denote by R{z) the unitary operator with the same kernel 
acting from L^(IR^) to M^^'^(IR^). We need two auxiliary spaces, "Ho £'^(1R) and 
Hi := C", and the corresponding trace maps tj : M^^'^(IR^) — Hj which act as 

Tqi/J := fs , Tllp := tn= (V" ■ • • > ) , 

respectively; they allow us to define the canonical embeddings of R(z) to Hi, i.e. 

R,l(z) = nR{z) ■.L^^n^, Ru{z) - [Il^L{z)]* -.n^^L^ 

and Rji(z) = Tj'RLi(z) : Hi — >■ Jij, all expressed naturally through the free Green's 
function in their kernels, with the variable range corresponding to a given Jii- The 
operator- valued matrix r(z) — [rij{z)] : 'Ho © Hi — > "Ho ® "Hi is defined by 



roo(z)/ 

Tii{z)ip 



= —Hij{z)g for i j and g EHj , 

RooW] / if /eHo, 



09 for C9 G 'Hi 



where s^, (z) — Pi + s{z) :— (3i + ^^(In ^ — and —'0(1) is the Euler number. 

For z from p(Ha^i3) the operator T{z) is boundedly invertible. In particular, roo(2;) 
is invertible which makes it possible to employ the Schur reduction procedure one 
more time and to define the map D{z) ■.'Hi ^'Hi by 

(9.2) D{z) = Tn{z) - Tw{z)Too{z)-'Toiiz) . 

We call it the reduced determinant of F; it allows us to write the inverse of F(z) as 
[F(z)]"i : Ha® Hi ^ Ho® Hi with the 'block elements' defined by 

[F(z)]7/ D{z)-\ 

[T{z)]-^ = Too{z)-' +roo{z)-'Toi{z)D{z)-^Tio{z)roo{z)-\ 

[r(z)]o7 = -Foo(z)-^Foi(z)i?(z)-i, 

[r(z)]7o' = -D{z)-'Tio{z)Tooiz)-' ; 

in the above formulae we use the notation Tij{z)^^ for the inverse of Tij{z) and 
[F(z)]~"'^ for the matrix element of [F(z)]^^. 

Before using this to express the resolvent Ra.p{z) := (Ha^is — z)^^ we introduce 
another notation which allow us to write Ra^piz) through a perturbation of the 
'line only' Hamiltonian Ha d escribing t he system without the point interactions, 
i.e. Pi — oo, i — 1, . . . ,n. By (BEKS94j the resolvent of Ha is equal to 

Raiz) = R{z) + Rlo(.z)Too'Rol{z) 



44 



PAVEL EXNER 



forz e C\[-ia2,cx)). We define the map Ra;Li(z) : Hi L'^{R^) hy Ila-iL{z)->P 
Ra{z)ip \n and Ka-iLiz) : L^(M^) — > "Hi as its adjoint, Rq;li(z) :— R*.;^^(z). 
The resolvent difference between Ha^p and Ha is then given by Krein's formula 
[AGHH] ■ A straightforward computation |EK04| yields now the following result. 

Theorem 9.1. For any z e pili^.fi) with Imz > we have 
1 

The obtained resolvent expressions allow us to investigate various spectral proper- 
ties of the operator Ha^p |KK04j; here we concentrate only on those related to the 
subject of the paper, namely to perturbations of embedded eigenvalues. 

9.2. Resonance poles. The mechanism governing resonance and decay phe- 
nomena in this model is the tunneling between the points and the line. This inter- 
action can be 'switched off' if the line is removed, in other words, put to infinite 
distance from the points. Consequently, the 'free' Hamiltonian Hp :— i?o,/3 has the 
point interactions only. It has m eigenvalues, 1 < m < n of which we assume 

(9.3) - < ei < ••• < e„ < 0, 

i.e. that the discrete spectrum of Hp is simple and contained in (the negative 
part of) a{Ha) — (Tac{Ha,p) — [—jo^joo); this can be always achieved by an 
appropriate choice of the configuration of the set 11 and the coupling parameters /?. 
Let us specify the interactions sites by their Cartesian coordinates, y*-*-* = (0^,0^). 
It is also useful to introduce the notations a = (ai, a„) and dij = ly*^*) — y^^^j for 
the distances between the point interactions. 

Resolvent poles will be found through zeros of the operator-valued function 



(9.2), more exactly, through the analytical continuation of D{-) to a subset ri__ of 
the lower halfplane across the segment (— ^a^, 0) of the real axis, in a similar way 



to what we did for Friedrichs model using formula (3.3). For the sake of definiteness 
we employ the notation D{-Y''\ where / — --1,0, 1 refers to the argument z from 
the segment {—ja'^,0), and the upper halfplane, Imz > 0, respectively. Using 
the resolvent formula of the previous section we see that the first component of the 
operator- valued function Z)(-)''Ms an n x n matrix with the elements 

rii;jfe(z)(') = -(1 - S^k) ^Ko{d^kV^) + Sjk{l3j + l/27r(ln^~ ^(1))) 
for all the /. To find an explicit form of the second component let us introduce 

_ ta (a-2^(z-t)l/^)e»(--*)'^'(l"'l+l°^l) ,ti/2(,._,.) 

and //^^(A,t) := lim,,_i.o+ /iij(A + iri,t). Using this notation we can rewrite the 
matrix elements of (FiorgQ^Foi)^''' (z) appearing in (9.2) in the following form. 



Jo t — A — 

00 



t — z — \a 



= W , \ „,2 + (? - l).9a,.j (^) for / = 1, 
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where V indicates again principal value of the integral, and 

„ . .C2^ — ^'^ -a(|a.| + |a,|)/2 i(z+oP / i)^'\c,^c,) 

^"''^^ ^ (z + «2/4)l/2*' 

Using these formulae one has to find zeros of det_D(-)'^~-'^'; we shall sketch the argu- 
ment referring to |EK04|, IEIK07| for details. We have mentioned that resonances 
are caused by tunneling between the parts of the interaction support, hence it is 
convenient to introduce the following reparametrization. 



6(a) = (61(a), 6„(a)) with h{a) :-- 



and to put 7'i{b{a),z) := det £'(z)'^~^^. Since the absence of the line-supported 
interaction can be regarded as putting the line to an infinite distance from the 
points, it corresponds to 6 = in which case we have 77(0,2) — detrii(2;) and the 
zeros are nothing else than the eigenvalues of the point-interaction Hamiltonian 
in other words, 77(0, e^) = , i = 1, . . . ,m. Then, in analogy with Sec. [sj we have 
to check that the hypotheses of the implicit-function theorem are satisfied which 
makes it possible to formulate the following conclusion. 

Proposition 9.2. The equation •q{h,z) = has for all the hi small enough 
exactly m zeros which admit the following weak-coupling asymptotic expansion, 

Zi{b) = e^ + 0{\b\)+iO{\b\) where |6| := max 6^ . 

l<2<m 

This result is not very strong, because it does provides just a bound on the 
asymptotic behavior and it does not guarantee that the interaction turns embedded 
eigenvalues of i/^ into true resonances. This can be checked in the case n — 1 
|EK04j but it may not be true already for ri = 2. The simplest example involves a 
pair of point interactions with the same coupling placed in a mirror-symmetric way 
with respect to E. The Hamiltonian can be then decomposed according to parity, 
its part acting on functions even with respect to E has a resonance, exponentially 
narrow in terms of the distance between the points and the line, while the odd one 
has a embedded eigenvalue independently of the distance. On the other hand, if 
the mirror symmetry is violated, be it by changing one of the point distances or 
one of the coupling constants, the latter turns into a resonance and one can derive 



a weak-perturbation expansion [EK04] in a way similar to those of Sec. 4.1 

Let us also note that the explicit form of the resolvent given in Theorem |9.1| 
makes it possible to find the on-shell S-matrix from energies from the interval 
{—jO^, 0), that is, for states traveling along the 'wire', and to show that their poles 
coincide with the resolvent poles; for 777 = 1 this is done in |EK04| . 

9.3. Decay of the 'dot' states. The present model gives us one more op- 
portunity to illustrate relations between resonances and time evolution of unstable 
systems, this time on bound states of the quantum 'dots' decaying due to tunneling 



between them and the 'wire'. By assumption (9.3) there is a nontrivial discrete 
spectrum of Hp embedded in (— |a^,0), and the respective eigenfunctions are 

m I 

M^) = E^i'Vi'^(^) , J = 1, • • • ,m , 4'\x) := J M^\x - y«|) , 
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where in accordance with |AGHHl Sec. II. 3] the coefficient vectors G C™ solve 
the equation rii(ej)c?*^-'' = and the normahzation condition = 1 gives 

m 2—1 

|d(^)p + 2Re5:^d(^)d^)(0p\0i^)) = l. 

'i=2 k=l 

In particular, d^^-* = lifn = TO = l;ifm>l and the distances between the points 
of n are large, the natural length scale being given by (— ej)^^/^, the cross terms 
are small and the vector lengths are close to one. 

Let us now identify the unstable system Hilbert space Hu = E^L'^{M.'^) with 
the span of the vectors ipi, . . . , ipm- The decay law of the system prepared at the 



initial instant t at a state ip € is according to (2.1 ) given by the formula 

We are particularly interested in the weak- coupling situation which in the present 
case means that the distance between S and 11 is a large at the scale given by 
(— Cm)"^''^. Let us denote by Ej the one-dimensional projection associated with 
the eigenfunction tpj, the one can make the following claim |EIK07| . 

Theorem 9.3. Suppose that Ha.p has no embedded eigenvalues. Then, with 
the notation introduced above, we have in the limit \b\ — ^ 0, i.e. dist (S,n) oo 

||£;,e-*^°-'^Vj-e-'"'Vjll ^0, 

pointwise in t £ (0, oo), which for n ~ 1 implies jP-^^ (t) — e^^™^i*| — > as |6| — > 0. 

Let us add a couple of remarks. The result implies more generally that for large 
values of dist (S, 11) the reduced evolution can be approximated by a semigroup. 
On the other hand, despite the approximately exponential decay in the case n = 1 
the lifetime defined as T^j — Xi -Pvi (^) "^^ diverges; the situation is similar to those 
mentioned is Sections 3] and |5.1[ the operator Ha^p has a bound state which is 
not exactly orthogonal to t/ji for b ^ 0, cf. |EK04j . hence \imt-^^ P^-^{t) ^ 0. 
Furthermore, the decay of the 'dot' states in this model offers a possibility to 
compare the 'stable' dynamics, i.e. evolution of vector in Hu governed e"'^''*, with 
the Zeno dynamics obtained from e"*^"'-^* by permanent observation, cf. |EIK07j 
for details. Finally, let us finally mention that a related model with a singular 
interaction in supported by a line and a circle and resonances coming from a 
symmetry violation has been investigated recently in jKol2| . 

10. Generalized graphs 

In the closing section we will mention another class of solvable models in which 
resonances can be studied, which may be regarded as another generalization of the 
quantum graphs discussed in Section [7j What they have in common is that the 
configuration space consists of parts connected together through point contacts. In 
the present case, however, we consider parts of different dimensions; for simplicity 
we limit ourselves to the simplest situation when the dimensions are one and two. 

10.1. Coupling different dimensions. To begin with we have to explain 



how such a coupling can be constructed. The technique is known since |ES87 



we demonstrate it on the simplest example in which a halfline lead is coupled to a 
plane. In this case the state Hilbert space is L^(M~) ® L^(IR^) and the Hamiltonian 
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acts on its elements (^'""'^ ) (belonging locally to W^'^) as (_^^'°'"' ); to make 
such an operator self-adjoint one has to impose suitable boundary conditions which 
couple the wave functions at the junction. 

The boundary values to enter such boundary condition are obvious on the 
lead side being the columns of the values ?/'iGad(0-l-) and '0i'ead(O+)- other 
hand, in the plane we have to use generalized ones analogous to those appearing 
in the first relation of (9.1 ). If we restrict two-dimensional Laplacian to functions 
vanishing at the origin and take an adjoint to such an operator, the functions in 
the corresponding domain will have a logarithmic singularity at the origin and the 
generalized boundary values will be the coefficients in the corresponding expansion, 

V'piano(a;) = io('/'piano) In |a;| + ii('(/'piano) + o{\x\) ; 
using them we can write the sought coupling conditions as 

^^p^^ V'l'ead(0+) = AV'lead(0+)+27rCLo(V'pla„e), 

il('0planc) = C'V'lead(O-H) Dio(V'planc) , 

where A, Z? e M and C is a complex number, or more generally 

V'lead(0+)\ , ^^CadlO- 



^ V^o(V'plane)/ ^ V^l (V'planc)y 



with appropriately chosen matrices A,B in analogy with (7.2), however, for our 
purpose here the generic conditions (10.1) are sufficient. 

As in the case of quantum graphs the choice of the coupling based on the 
probability current conservations leaves many possibilities open and the question 
is which ones are physically plausible. This is in general a difficult problem. A 
natural strategy would be to consider leads of finite girth coupled to a surface and 
the limit when the transverse size tends to zero. While for quantum graphs such 
limits are reasonably well understood nowadays |Gr08l IEP091 ICETlOj , for mixed 
dimens ions the current knowledge is limited to heuristic results such as the one in 
[ES97j which suggests that an appropriate parameter choice in (10.1) might be 

(10.2) A^^, B=.l^, C--^, D = -\np, 

where p is the contact radius. At the same time, other possibilities have been con- 
sidered such as the simplest choice keeping just the coupling term, ^ = Z? = 0, or an 
indirect approach based on fixing the singularity of the Hamiltonian Green's func- 
tion at the junction which avoids using the coupling conditions explicitly |Ki97j . 

While the example concerned a particular case, the obtained coupling condi- 
tions are of a local character and can be employed whenever we couple a one- 
dimensional lead to a locally smooth surface. In this way one can treat a wide class 
of such systems, in particular to formulate the scattering theory on configuration 
spaces consisting of a finite numbers of manifolds, finite and infinite edges — one 
sometimes speaks about 'hedgehog manifolds' — cf. |BG03| . 

Before turning to an example of resonances on such a 'manifold' let us mention 
that while the system of a plane and a halfline lead considered above has at most 
two resonances coming from the coupling, one can produce an infinite series of them 
if the motion in the plane is under influence of a magnetic field. The same is true 
even if Laplacian is replaced by a more complicated Hamiltonian describing other 
physical effects such as spin-orbit interaction — cf . [CEllj . 
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10.2. Transport through a geometric scatterer. Let us look in more de- 
tail into an important example in which we have a 'geometric scatterer' consisting 
of a compact and connected manifold fl, which may or may not have a boundary, 
to which two semi-infinite leads are attached at two different points xi, X2 from the 
interior of ^l. One may regard such a system as a motion on the line which is cut 
and the loose ends are attached to a black-box object which can be characterized 
by the appropriate transfer matrix, ("/(g^^)) = i( "/(q""*)) • To find the latter one 
has to fix the dynamics: we suppose that the motion on the line is free being de- 
scribed by the negative Laplacian, while the manifold part of the Hamiltonian is 
Laplace-Beltrami operator on the state Hilbert space iF'iQ) of the scatterer; they 
are coupled by conditions ( 10.1 1 with the coefficients indexed by j — 1,2 referring 
to the 'left' and 'right' lead, respectively. 

We need the Green function G(., .; fc) of the Laplace-Beltrami operator which 
exists whenever the does not belong to the spectrum. Its actual form depends on 
the geometry of fl but the diagonal singularity does not: the manifold n admits in 
the vicinity of any point a local Cartesian chart and the Green's function behaves 
with respect to those variables as that of Laplacian in the plane, 

G(x, y; fc) = - ^ In \x-y\ + 0(1) , \x~y\ ^ . 

Looking for transient solutions to the Schrodinger equation at energy k^, we note 
that its manifold part can be written as u{x) = aiG{x,xi;k) + a2G{x,X2;k), 
cf. |Ki97j . which allows us to find the generalized boundary values 

Loixj) = , Li{xj) = aj^{xj,k) + a3^jG{xi, X2. k) 

ZTT 

for j — 1, 2, where we have employed the regularized Green's function at Xj, 

In \x—x- 



(10.3) ^{xj;k) := lim 



G{x, Xj] k) -I- 



27r 



Let Uj be the wave function on the j-th lead; using the abbreviations Uj, u'^ for 
its boundary values we get from the conditions (10.11 a linear system which can 
be easily solved [ETVOl] : it yields the transfer matrix in terms of the quantities 
Zj := ^ + and A := — Z1Z2 , where £,j := ^{xj;k) and g := G{xi,X2',k). 
The expression simplifies if the couplings are the same at the two junctions; then 
det L — 1 and the transfer matrix is given by 

Z2 + ~2j^ \ 




- AiZi + Z2)- ^A + Zi J 

From here one can further derive the on-shell scattering matrix [ETVOl] . in par- 
ticular, the reflection and transmission amplitudes are 

L21 + ik{L22-Lii) + k^Li2 ^ 2ik 

L21 - ik{L22+Lii) - k^Li2 ' L21 - ik{L22 + Lii) - k'^Li2 ' 

they naturally depend on k through ^ and 5, and satisfy |rp-f |tp — 1. To find these 
quantities for a particular VI we may use the fact that it is compact by assump- 
tion, hence the Laplace-Beltrami operator on it has a purely discrete spectrum. We 
employ the eigenvalues, {\n\'^=i , numbered in ascending order and with the multi- 
plicity taken into account, corresponding to eigenfunctions {4>n}'^=i which form an 
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orthonormal basis in L (fi). The common Green's function expression then gives 



g{k) - ^ 



^ A„- fc2 



while the regularized value ( |lQ.3[ ) can be expressed |ES97 as 



=1 



where the series is absolutely convergent and the constant c(0) depends on the 
manifold G. Note that a nonzero value of c(n) amounts in fact just to a coupling 
parameter renormalization: Dj has to be changed to Dj + 27rc(f7) . 

Several examples of such a scattering has been worked out in the literature, 



mostly for the case when fi is a sphere. If the coupling is chosen according to ( 10.2 ) 
and the leads are attached at opposite poles, the transmission probability has reso- 
nance peaks around the values A„ where the transmission probability is close to one, 
and a background, dominating at high energies, which behaves as 0{k~'^{\nk)~^), 
cf. [ETVOlj . Similar behavior can be demonstrated for other couplings at the junc- 
tions |Ki97| : the background suppression is faster if the junctions are not antipolar 
|BGMP02] . Recall also that this resonance behavior is manifested in conductance 
properties of such systems as a function of the electrochemical potential given by 
the Landauer-Biittiker formula, see e.g. |BGMP02j . 

10.3. Equivalence of the resonance notions. Let us return finally to a 
more general situatiorj^ and consider a 'hedgehog' consisting of a two-dimensional 
Riemannian manifold Q, compact, connected, and for simplicity supposed to be 
embedded into M^, endowed with a metric grs, to which a finite number Uj of 
halfline leads is attached at points Xj, j = 1, . . . , n belonging to a finite subset {xj} 
of the interior of fl; by M — J^j "-j '^6 denote the total number of the leads. The 

Hilbert space will be correspondingly H = L'^{n, v/|g[da;) ® 0fl^ L'^{Rf). 

Let Hq be the closure of the Laplace-Beltrami operator —g^^^^drig^^^g^^dg) 
defined on functions from C^{il); if dil 7^ we require that they satisfy there 
appropriate boundary conditions, either Neumann/Robin, (9„ -1- j)f\dn = 0, or 
Dirichlet, /jgo = 0. The restriction Hq of Hq to the domain {/ e D{Ho) : f{xj) = 
0, J = 1, . . . , n} is a symmetric operator with deficiency indices (n, n). Further- 
more, we denote by Hi the negative Laplacian on L2(IR^^) referring to the i-th 
lead and by H'^ its restriction to functions which vanish together with their first 
derivative at the halfline endpoint. The direct sum H' — Hq ® ® • • • ® H'j^j is 
obviously a symmetric operator with deficiency indices {n + M, n + M). 

As before admissible Hamiltonians are identified with self-adjoint extensions 



of the operator H' being described by the conditions (7.2 1 where U is now an 
(n+M) X (n+M) unitary matrix, / the corresponding unit matrix, and furthermore, 
* = (■^ia(/),-----^i,n(/),/i(0),...,/n(0))'^ and is the analogous column of 
(generalized) boundary values with Lij{f) replaced by Loj{f) and /i(0) by fj{0), 
respectively. The first n entries correspond to the manifold part being equal to the 



Considerations of this section follow the paper IEL12I . Similarly one can treat 'hedgehogs' 
with three-dimensional manifolds, just replacing logarithmic singularities by polar ones. 
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appropriate coefficients in the expansion of functions / G D(Hq) the asymptotic 
expansion near xj, namely f{x) — LQj{f)FQ{x, Xj) + Lij{f) + 0{r{x, Xj)), where 

Fo{x,Xj) = — lnr{x,Xj) 

Zn 

with r{x, Xj) being the geodetic distance on fl; according to Lemma 4 in [BG03] 52 
is a continuous functions of x with qi{xj, Xj) = 1. The extension described by such 
conditions will be denoted as Hi/. We are naturally interested in local couplings; 
in analogy with considerations of Section [7.1| we can work with one 'large' matrix 
U and encode the junction geometry in its block structure. 

Another useful thing we can adopt from the previous discussion is the possibility 
to employ conditions (Uj{k) — I)dj(f) + i{Uj{k) + I)cj{f) = on the manifold 
itself with the effective, energy-dependent coupling described by the matrix 

Uj{k) = Uij - (1 - fc)L/2j[(l - k)U4j -{k + l)ir^U3j 

at the j-th lead endpoint, where Uij denotes top- left entry of Uj, U2j the rest of 
the first row, U^j the rest of the first column and U^j is rij x Uj part corresponding 
to the coupling between the leads attached to the manifold at the same point. 

To find the on-shell scattering matrix at energy fc^ one has to couple solutions 
aj{k)e~'^^^ -\- hj{k)e^^^ on the leads to solution on the manifold and to look at 
the continuation of the result to the complex plane. On the other hand to find 
the resolvent singularities, we can again employ the complex scaling and to find 
complex eigenvalues of the resulting non-self adjoint operator. In both cases we 
need the solution on the manifold; modifying the conclusions of |Ki97| mentioned 
above we can infer that it has to be of the form f(x, k) = CjG{x, Xj] k), 

Consider first the scattering resonances. Denoting the coefficient vector of 
f{x,k) as c and using similar abbreviations a for the vector of the amplitudes of 
the incoming waves, (ai(fc), . . . , aM{k))'^ , and b for the vector of the amplitudes of 
the outgoing waves, one obtains in general a system of equations, 

A{k)a + B{k)h + C{k)c = 0, 

in which A and B are (n+M) x M matrices and C is (n+M) x n matrix the elements 
of which are exponentials and Green's functions, regularized if necessary; what is 
important that all the entries of the mentioned matrices allow for an analytical 
continuation which makes it possible to ask for complex k for which the above 
system is solvable. For /cq ^ M the columns of C(fco) are linearly independent and 
one can eliminate c and rewrite the above system as 

A{ko)a + B{ko)h = 0, 

where A{ko) and B{ko) are M x M matrices the entries of which are rational 
functions of the entries of the previous ones. If detA{ko) = there is a solution 
with b = 0, and consequently, fcg should be an eigenvalue of H since Imfco < 
and the corresponding eigenfunction belongs to L^, however, this contradicts to the 
self-adjointness of H. Next we notice that the S-matrix analytically continued to 
the point ko equals —B{ko)~^ A{ko) hence its singularities must solve det B{k) = 0. 

On the other hand, for resolvent resonances we use exterior complex scaling 
with arg6' > argfco, then the solution aj{k)e~^^^ on the j-ih lead, analytically 
continued to the point k = ko, is after the transformation by U0 exponentially 
increasing, while bj{k)e^^^ becomes square integrable. This means that solving in 
the eigenvalue problem for the complex-scaled operator one has to find solutions 



SOLVABLE MODELS OF RESONANCES AND DECAYS 



51 



of the above system with a = which leads again to the condition deti3(A;) = 0. 
This allows us to make the following conclusion. 

Theorem 10.1. In the described setting, the hedgehog system has a scattering 

resonance at with Im fcg < and /cq ^ M iff there is a resolvent resonance at fco • 
Algebraic multiplicities of the resonances defined in both ways coincide. 
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